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Abstract 

The first flavor-tagged amplitude analysis of the decay D° to the self-conjugate final state 
K + K~it + 'k~ is presented. Approximately 3000 signal decays are selected from data acquired 
by the CLEO II.V, CLEO III, and CLEO-c detectors. The three most significant amplitudes, 
which contribute to the model that best fits the data, are (f>p°, -fCi(1270) ± -fC =F , and non-resonant 
K + K~7r + 7r~ . Separate amplitude analyses of D° and D° candidates indicate no CP violation 
among the amplitudes at the level of 5% to 30% depending on the mode. In addition, the sen- 
sitivity to the CP-violating parameter 7/(^3 from a sample of 2000 B + — > D°(K + K~ir + ir~)K + 
decays, where D is a D° or D°, collected at LHCb or a future flavor facility, is estimated to be 
(11.3 ± 0.3)° using the favored model. 



*Present address: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 



A rich variety of interesting physics may be explored by investigating the decay D° — > 
K~ K + 'k + 'k~ . Study of the relative contribution of the intermediate resonances participating 
in the decay can help in the understanding of the behavior of the strong interaction at low 
energies. The mode is also of interest for its application in CP violation studies, both for 
improving the knowledge of the CKM unitarity triangle, and in probing for new physics 
effects through direct CP violation searches in D meson decays. Although results on the 
resonant structure of the decay have been reported by the E791 [l[ and FOCUS [2[ collab- 
orations, studies of higher precision are required. In particular, neither of these previous 
analyses differentiated between D° and D° decays. 

An important goal in flavor physics is the precise determination of the CKM unitarity 
triangle angle 7 (denoted by others as ^3), the phase of V c b relative to V u b- This parameter 
can be measured through the study of interference effects in the decay B^ — > D°K ± . Here 
D° indicates either a D° or a D° meson decaying to a hadronic final state. Experimentally, 
in order to obtain the best possible knowledge of 7 it is important to make use of as many 
D° decay modes as possible. The decay D° — > K + K~7f + 7t~ has been noted as being of 
potential interest in this respect [jjj], especially since large numbers of B ± — > D°K ± decays 
involving this mode will be collected by the LHCb collaboration. The existing knowledge of 
the substructure of D° — > K + K^-k + -k^ is however inadequate to make a reliable assessment 
of the potential sensitivity to 7, and so improved information is required. 

A search for direct CP violation in singly- Cabibbo-suppressed (SCS) charm decays is a 
promising method to test for the contribution of new physics, which in several plausible 
scenarios could lead to (9(1%) effects, an order of magnitude higher than is expected in the 
Standard Model Q. Evidence of CP violation has recently been reported in two-body SCS 
decays j5[, hence it is important to look elsewhere. The decay D° — > K + K~tt + tt~ , with 
a rich structure of intermediate resonances, is a suitable mode in which to perform such 
a search. CP violation studies involving D° — > K + K~7r + 7r~ have been conducted by the 
FOCUS [6[ and BABAR j3] collaborations, using the method of T-odd correlations, and 
null results have been reported. A CP violation search made with an amplitude analysis 
remains valuable however, as it probes each intermediate resonance of the decay separately, 
and hence can expose effects which may be diluted or concealed by the more inclusive T-odd 
correlation approach. 

This paper describes a flavor-tagged amplitude analysis of the decay D° — > K + K^-k + -k^ 
made using data collected by the CLEO collaboration in e + e~ collisions at the Cornell 
Electron Storage Ring (CESR). An amplitude model is constructed and a CP violation 
study performed. The model is also used to assess the potential sensitivity of the decay in a 
future — > D°K ± 7 measurement. The data analysed were collected with several different 
configurations of the CLEO detector, and at different center-of-mass energies. Comparison of 
the results obtained for each data set provides a powerful test of systematic robustness. One 
sample consists of CP-tagged decays from CLEO-c running at the ^(3770). These events 
provide unique access to the strong-phase differences between the intermediate resonances 
in a manner which was not available to previous studies. 

The paper is organized as follows. Section |TT] discusses the data sets used in the analysis. 
Section 0TI] describes the amplitude fit procedure and the development of the resonance 
model. Section IIVI presents the final model, summarizes the systematic uncertainties and 
gives the result of the CP violation test. The sensitivity of the decay D° — > K + K~7r + 7r~ in 
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a measurement of 7 with B — > D°K ± decays is considered in Sec. |V] and conclusions are 
given in Sec. ED 

II. DATA SET AND EVENT SELECTION 

The data analysed in this paper were produced in symmetric e + e~ collisions at CESR be- 
tween 1995 and 2008, and collected with three different configurations of the CLEO detector: 
CLEO II.V, CLEO III, and CLEO-c. 

In CLEO II.V [§[ tracking was provided by a three-layer double-sided silicon vertex detec- 
tor, and two drift chambers. Charged particle identification came from dE/dx information 
in the drift chambers, and time-of-flight (TOF) counters inserted before the calorimeter. 
For CLEO III [9[ a new silicon vertex detector was installed, and a ring ima ging Cherenkov 
(RICH) detector was deployed to enhance the particle identification abilities 1 10] . In CLEO- 



c, the vertex detector was replaced with a low- mass wire drift chamber [111]. A supercon- 
ducting solenoid supplied a 1.5 T magnetic field for CLEO II.V and III, and 1 T for CLEO-c 
operation, where the average particle momentum was lower. In all detector configurations 
neutral pion and photon identification was provided by a 7800-crystal Csl electromagnetic 
calorimeter. 

Four distinct data sets are analysed in the present study: 

1. approximately 9 fb _1 accumulated at \fs ~ 10 GeV by the CLEO II.V detector; 

2. a total of 15.3 fb~ x accumulated by the CLEO III detector in an energy range \fs = 
7.0 - 11.2 GeV, with over 90% of this sample taken at ^/s = 9.5 - 10.6 GeV; 

3. 818 pb" 1 collected at the ^(3770) resonance by the CLEO-c detector; 

4. a further 600 pb" 1 taken by CLEO-c at -y/s = 4170 MeV. 

These samples are referred to as the CLEO II.V, CLEO III, CLEO-c 3770 and CLEO-c 

4170 data sets, respectively. 

The analysis considers two classes of signal decays, for both of which information on the 
quantum numbers of the meson decaying to the signal mode is provided by an event tag. 

• Flavor-tagged decays are selected from the CLEO II.V and CLEO III data sets, in 
which the flavor of the decaying meson is determined by the charge of the 'slow pion', 
7r s , in the D* + — > P°7r+ decay chain. Flavor-tagged decays are also selected from the 
two CLEO-c data sets, where here the tag is obtained through the charge of a kaon 
associated with the decay of the other D meson in the event. 

• CP-tagged decays are selected in the CLEO-c 3770 data set alone. In ^(3770) 
decays the D — D pair is produced coherently. Therefore, the CP of the signal D can 
be determined if the other D meson is reconstructed in a decay to a CP-eigenstate. 
Useful information is also obtained if the tagging meson is reconstructed decaying into 
the modes Kgii + TT^ or A"£7r + 7r _ , for which the relative contribution of CP-even and 
CP-odd states is known [12j. 

Detector response is studied with GEANT-based 13[ Monte Carlo simulations of each 
detector configuration, in which the Monte Carlo events are processed with the same recon- 
struction code as used for data. 
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A. Flavor-tagged CLEO II. V and CLEO III samples 



Selections are run on both the CLEO II. V and CLEO III data sets to identify events 
containing the fully reconstructed decay chain D* + — > -D°7r+, D° — > K + K~7t + tt~ . These 
selections are not identical for the two data sets on account of the different detector responses. 

Requirements are first placed on the attributes of the individual charged tracks used in the 
reconstruction. The tracks must be well-measured and satisfy criteria based on fit quality. 
They must also be consistent with coming from the interaction point in three dimensions. 
In the CLEO III analysis, the polar angle of each considered track is required to satisfy 
| cos#| < 0.9. The momentum of the slow pion candidate must be above 100 (150) MeV/c 
in the CLEO II.V (CLEO III) analysis and below 500 MeV/c, and that of the other final 
state particles must be between 200 MeV/c and 5000 MeV/c. 

Particle identification information plays an important role in the selection. In the 
CLEO III analysis, candidate tracks with momentum above 500 MeV/c are classified as 
kaons if they have at least three associated photons in the RICH detector, and a ring fit 
to the photon hits indicates that the kaon hypothesis is more probable than that of the 
pion. Lower momentum tracks, and those tracks lying outside the angular acceptance of the 
RICH, are identified as kaons if they have a dE/dx value within 2.1 o of that expected for a 
true kaon. Pion candidates are required to have a dE/dx value within 3.2 a of that expected 
for a true pion. In the CLEO II.V study, tracks with dE/dx information are identified as 
kaon candidates if they have a dE/dx value lying within 2.1 (2.5) a of that expected for a 
true kaon (pion). When TOF information is available, both kaon and pion candidates are 
required to lie within 2.5 a of their expected value. 

A possible background to the signal arises from the decay D° — > K s (-k + n~)K + K~ , and 
so a K s veto procedure is performed. If the two pions have an invariant mass compatible 
with a Kg decay the event is rejected if the flight distance of the K s candidate from the 
interaction point, normalized by the assigned error, is greater than two, or either pion has 
an impact parameter in the transverse plane greater than 1.5 mm. 

After all selection criteria 3.0% of the remaining events in the CLEO II.V sample are 
found to have more than a unique pair of D* + and D° candidates. If there is more than 
one D* + candidate, the D* + chosen chosen is the one which has a D* + — D° invariant mass 
difference closest to the expected value. If there is more than one D° candidate, the one 
that is chosen is the D° for which the dE/dx information of the four daughter tracks best 
matches the signal hypothesis. In the CLEO III sample only 0.3% of events contain multiple 
candidates; these events are discarded. 

Two kinematic fits are performed to the decay with the constraint that the four tracks 
from the D° meson candidate originate from a common vertex, and also that the D° candi- 
date and the slow pion from the D* + originate from another common vertex. Loose criteria 
are placed on the quality of these fits. The analysis is optimized to favor charm mesons pro- 
duced in the primary interaction, rather than B meson decay, and so it is required that the 
D* + momentum is at least half the maximum kinematically allowed value. This requirement 
suppresses combinatoric background. 

Figure [1] shows the spectrum of the D° candidate invariant mass, mo, and that of the 
invariant mass difference, Am, between the D* + and D° candidates for the CLEO II.V 
selection after the vertex-constrained fit. In Fig. [2] are shown the equivalent plots for the 
CLEO III analysis. 

The rri£, and Am distributions are fit with single (sum of bifurcated) Gaussians for the 
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FIG. 1: (a) Am (after mu selection) and (b) ran (after Am selection) distributions of CLEO II. V 
D° candidates. 




FIG. 2: (a) Am (after mp selection) and (b) mo (after Am selection) distributions of CLEO III 
D° candidates. 



signal peaks in the CLEO II. V (CLEO III) sample, and empirical functions to describe the 
background. A signal region is defined as being within ±5.0 MeV/c 2 (^i! 3 2 MeV/c 2 ) of the 
world average value of the D° mass [3] for mo in the CLEO II. V (CLEO III) case, and within 
±0.80 MeV/c 2 (l^so MeV/c 2 ) of the world average value of the D* + — D° mass difference for 
Am. An analysis of the 7r + 7r~ invariant mass spectrum, detailed in Sec. IHH indicates that 
the peaking background from D° — > K®K + K~ events in the sample is negligible, and so the 
overall signal yield and background level can be extracted directly from these fits. For both 
data sets the fits to the and Am spectra give consistent results. In total 279 events 
are selected in the signal region from the CLEO II. V data set, with an estimated purity of 
(74 ±3)%. In the CLEO III analysis 1225 events are selected with a purity of (89.2 ±0.4)%. 

In order to learn about the characteristics of the contamination in the signal region, 
dedicated background samples are also selected. These are taken from three separate regions 
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of uid — Am space: one with the same acceptance on mo as the signal region, but with 
148.5 < Am < 155.5 MeV/c 2 , and two others with the same acceptance on Am as for the 
signal sample, but with 1.7 < mjj < 1.8 GeV/c 2 and 1.9 < mo < 1.95 GeV/c 2 respectively. 
Simulation studies indicate that these regions contain negligible contributions from true 
signal decays, and that the attributes of the selected events are representative of those of 
the background events in the signal sample. 

In order to improve the resolution of the four-momenta used in the amplitude analysis, 
the selected events in the signal sample are subjected to a kinematical fit in which the D° 
candidate daughter particles are constrained to originate from a common vertex, and to 
have an invariant mass equal to the nominal D° mass. The same procedure is applied to 
the selected events in the background sample. 

The performance of the flavor tag is calibrated in data using D° — > K~7t + 7t~tt + decays. 
These decays, accompanied by a slow pion, are selected with the same procedure as for the 
signal sample. By comparing the charge of the tagging slow pion with that of the slow 
pion from the fully reconstructed D decay it is possible to determine directly the mistag 
rate, that is, the fraction of occasions on which the tagging decision is incorrect. Small 
corrections are applied to account for the contribution from doubly-Cabibbo-suppressed 
decays. Simulation studies are used to validate that the mistag rate as determined by this 
procedure is consistent with that of D° — > K + K~tt + it~ decays. This study is only performed 
on the CLEO III sample, yielding a mistag rate of (0.64±0.05). In the subsequent amplitude 
analysis the same value is taken to apply for the CLEO II. V sample, and a relative 50% 
uncertainty assigned to account for this assumption. 



B. Flavor-tagged CLEO-c samples 

Flavor tagging is performed in both CLEO-c data sets by searching for another charged 
kaon in the event, in addition to those used in reconstructing the signal D decay. Such a 
tagging kaon originates from the decay of the other D meson. If this decay is assumed to 
be Cabibbo-favored, and there are no other additional kaons in the event, then the charge 
of the tagging kaon indicates the flavor of the decaying meson, and hence the flavor of the 
signal decay can be inferred. 



Standard CLEO-c selection criteria, as described in Ref. [15], are imposed on the tracks 
used in the D° reconstruction and for the tagging kaon. Events are only considered with 
a single tagging kaon candidate, on which a momentum cut is applied. Selecting higher 
momentum kaons is found to be advantageous both in enhancing the purity of the sample, 
and in suppressing events where the tagging decision is incorrect. The momentum of the 
tagging kaon is required to exceed 400 MeV/c in the CLEO-c 3770 data set, and 600 MeV/c 
in the CLEO-c 4170 data set. 

It is necessary to apply a more stringent K s veto to suppress D° — > KgK + K~ con- 
tamination than in the CLEO II. V and CLEO III selections. This is because at CLEO-c 
the D mesons are produced at or close to threshold, and hence the flight distance of any 
resulting Kg is lower. Therefore, events are rejected in which a K s candidate has a flight 
distance, normalized by the assigned uncertainty, of greater than one, or in which either of 
the daughter pions has an impact parameter in the transverse plane greater than 1 mm. 
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1. CLEO-c 3770 sample 



Two kinematical variables are defined: the beam-constrained candidate mass, 



m bc 



Pl/c 2 



where p_o is the momentum of the signal D candidate, and AE = E D — y/s/2, where Ed is 
the sum of the energies of the daughter particles of the signal D candidate. The distributions 
of mbc and AE are shown in Fig. [3]for kaon-tagged candidates in the CLEO-c 3770 data set. 
The signal decays peak at the nominal D° mass in m& c and zero in AE. In making the final 
selection a window of ±5 MeV/c 2 and ±15 MeV is placed around these expected values for 
rribc and AE respectively. A sample of 1396 events is selected in the signal region, of which 
14 contain two candidates. 

In this latter class of event only one candidate, chosen at random, is retained for sub- 
sequent analysis. Interpolating the results of fits to the sideband regions into the signal 
window indicates that the contamination from non-peaking background is at the level of 
(13.5 ± 0.5)%. The residual contamination from D° — > K$K + K~ decays is found to consti- 
tute (2.4 ± 0.4)% of the sample, as determined from the amplitude fit studies described in 
Sec. [fill 

In addition a sample of 763 events is selected for non-peaking background studies from the 
regions defined by -5 < (m bc - 1865MeV/c 2 ) < 5 MeV/c 2 and \AE ± 45MeV| < 30 MeV. 
A further sample of 445 events which fail the Kg veto, but pass all other signal selection 
criteria, are selected in order to characterize the residual -D ( 
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contammation. 
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FIG. 3: (a) m& c and (b) AE distributions of events passing the kaon-tagged D° — > K + K~it + it~ 
selection in the CLEO-c 3770 data set. Each distribution is plotted after applying the selection 
cut on the other variable. 



The performance of the flavor tag is calibrated in data using D° — > K~n + n~n + decays. 
These decays, accompanied by a tagging kaon, are selected with the same procedure as for 
the signal sample. The method is validated using simulated data and corrections are made 
for the doubly-Cabibbo-suppressed decays in the sample as for the CLEO III calibration. It 
is concluded that the mistag rate of signal events in data is (4.5 ± 0.5)%. 
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All selected events are subjected to a kinematical fit with the invariant mass of the 
candidate constrained to that of the D°, in order to provide the best possible resolution for 
the amplitude study. 

2. CLEO-c 4170 sample 

At \/s = 4170 MeV, pairs of charm mesons can be produced in a variety of configurations, 
including DD, D*D, D*D*, D*Dir, D+D~ and D+*D~. Several of these configurations may 
result in events which contain a D° accompanied by a D°, or a D° and afl". Depending on 
the production process and subsequent strong or electromagnetic decay, there will be one 
or more prompt pions or photons also present in the event. Even without reconstructing 
these additional particles it is possible to separate statistically the different production and 
decay categories, as they exhibit different distributions in m& c — AE space. This property 
has been exploited in Ref. [l6| to study charm production at these energies. 

D*D* events have the highest rate and intrinsic purity, and so these are isolated for 
the amplitude analysis. A variable AE S i g = a — b m bc is defined, where the coefficients 
a = 2.112 GeV and b = 1.12 c 2 are obtained from a fit to the distribution of simulated 
signal candidates, and events are selected in the region 2.005 < m bc < 2.040 GeV/c 2 and 
\AE — AE sig \ < 10 MeV. In addition, to suppress background further, a restriction is placed 
on the momentum of the D° candidate that it be above 450 MeV/c. Figure H] (a) presents the 
distribution of m bc with the cut on AE — AE sig applied, and Fig. H](b) the corresponding plot 
for AE — AE sig , after selecting on m^. A total of 739 events is selected, of which 5 contain a 
second candidate. In the case of these multiple candidate events, only one candidate, chosen 
at random, is retained for the subsequent analysis. 
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FIG. 4: Distributions of (a) m\, c (after AE — AE s ^ g selection) and (b) AE — A.E s i g (after m bc 
selection) for events passing the kaon-tagged D° — > K + K~ir + ir~ selection in the CLEO-c 4170 
data set. 

Two regions are selected to provide non-peaking background samples. These are defined 
by the criteria 2.005 < m bc < 2.040 GeV/c 2 , \(AE - AE sig ) ± 45MeV| < 20 MeV and 
Pd > 450 MeV/c. This selection yields a total of 978 background candidates. Simulation 
indicates that the density and nature of the candidates in this sample is compatible with 
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TABLE I: CP-tagged D -)■ K+R-tt+tt- events selected from the CLEO-c 3770 MeV data set. 



CP-even ta, 


gs 


C-r-oad taj 


?s 


Admixture 


ivlocle 


i ielu 


Moae 


i ielu 


Mode Yield 


K+R- 


11 




15 


K^TT+TT- 63 




8 


K S U{-K+TT-TT Q ) 


12 


^vr+TT- 87 


K° L 7r° 


9 


K° s cj)(K + K-) 


1 






7 




2 






14 


^^(vr+vr-vr ) 


1 










1 




Total 


49 




32 


150 



that of the contamination inside the signal window. It is thus estimated that the fraction of 
non-peaking background in the signal sample is (33 ± 1)%. The fraction of D° — > K®K + K~ 
decays in the signal sample is estimated to be (1.3 ± 0.4)%. 

The performance of the kaon flavor tag is determined using D° — > K~tt + 7t~tt + events 
in the same manner as for the 3770 MeV sample. The mistag rate is measured to be 
(7.5 ±0.7)%. 

As in the case of the 3770 MeV analysis, the selected candidates are refit with the mass 
of the D° imposed as a constraint. 

C. CP-tagged CLEO-c sample 

CP-tagged events are selected from the 3770 MeV CLEO-c data set in which both D- 
meson decays are reconstructed, one through its decay to K + K~tt + tt~ and the other to a 
CP-eigenstate. The latter decay provides a tag through which the CP eigenvalue of the 
signal decay can be determined. The signal decay is selected in the same manner as for the 
flavor-tagged sample, and is required to lie within ±5 MeV/c 2 of the nominal D° mass in 
rribc and within ±20 MeV of zero in AE. The selection criteria for the CP-tags are identical 
to those used in Ref. [17| . The CP-tags include the modes D — > K^tt and K^u, where the 
presence of the K\ meson is inferred from a missing- mass technique, having reconstructed all 
the other particles in the event. The number of selected candidates is presented in Table [B 

Events are also selected in which the tag is provided by either of the decays K° s ti + ti^ 
or K°7r + 7r~, reconstructed with the same requirements as in Ref. (l2| . These modes can 
be considered admixtures of CP-odd and CP-even eigenstates, and are exploited in the 
analysis thanks to available measurements of the strong-phase variation across the Dalitz 
space of each decay Table [J reports the yields for this category of event. 

The level of contamination is estimated from study of sidebands in the two-dimensional 
rribc space of the reconstructed D mesons, in the missing-mass sideband for the events con- 
taining K\ tag and using Monte Carlo to determine the contribution of peaking backgrounds. 
The fraction of contamination is found to be (35.9±5.4)% for the CP-even tags, (17.2±3.4)% 
for the CP-odd tags, (16.5 ± 2.3)% for the K^ir + ir~ tags and (19.0 ± 2.2)% for the K L n + 7r- 
tags. The sidebands do not provide a data set of sufficient size to allow the composition of 
the background to be studied. For this purpose, a simulated background sample is prepared. 

Multiple candidates occur in 1.3% of selected signal events. In these events only a single 



10 



candidate is propagated for further analysis, taking the signal-tag combination in which the 
sum of mbc values is closest to twice the nominal D° mass. In events where the tag involves 
a K° L , the candidate is selected for which m\, c of the signal decay is closest to the D° mass. 
All selected candidates are refit with a D° mass constraint. 

D. Summary of signal samples 

A summary of the tagged samples is presented in Table DTI A total of 3639 flavor-tagged 
events and 231 CP-tagged events are selected for the amplitude analysis studies, of which 
2959 and 181 events, respectively, are estimated to be true signal decays. 



TABLE II: Summary of signal samples. 



Sample 


Yield 


Purity [%] 


Mistag rate [%] 


Flavor tags 








CLEO II.V 


279 


74 ±3 


0.64 ±0.32 


CLEO III 


1225 


89.2 ±0.4 


0.64 ±0.05 


CLEO-c 3770 


1396 


84.1 ±0.6 


4.5 ±0.5 


CLEO-c 4170 


739 


65.7 ± 1.1 


7.5 ±0.7 


CP-tags 








CP-eigenstates 


81 


71.9 ±3.1 




Admixture 


150 


82.1 ± 1.5 





III. AMPLITUDE ANALYSIS 

The relative magnitudes and phases of the intermediate resonances are determined by a 
maximum-likelihood fit to the data selected. The formalism and composition of the like- 
lihood is described in Sec. IIII Al Sections IIIIBI and II II CI describe the calculation of the 
likelihood normalization and the goodness-of-fit, respectively. 

A. Likelihood fit 

The likelihood contains probability density functions (PDFs) for the signal amplitude 
model and background components as a function of particle four-momenta. These PDFs are 
modified to account for the variations in acceptance over phase space. The method to con- 
struct the PDFs is similar to that used to describe a decay to three final-state pseudoscalar 
mesons (three-body decay) performed previously by the CLEO Collaboration [l8j]. How- 
ever, a decay to four final-state pseudoscalar mesons (four-body decay) is more complicated 
due to the non-uniform phase space and the possibility of having two separate intermediate 
resonances contributing to the amplitude. 

The PDFs are functions of the four- momentum pj, of the D° and its decay products, 
where the index j =0, 1, 2, 3, and 4, corresponds to D°, K + , K~, 7r + , and 7r~ mesons, 
respectively. The four-body phase space function, R^(pj), is not uniform in any set of five 
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independent kinematic variables that are required to describe the phase space, unlike the 
uniform three-body phase-space function over the Dalitz plot. Therefore, an analytic form 
of Ri{pj) is used that accounts for the kinematic constraints among the four- momenta that 
arise from the invariant masses of the D° and its decay products [191 ] . 

There are several different amplitude types that can arise in four-body decays. 

• A single-resonance amplitude such as D° -»■ K*(892) K~'K + , K*(892)° -»■ K + -k~. 
Unless the resonance is a scalar, there are several possible amplitudes with the same 
intermediate resonance but with differing orbital angular momentum among the final 
state particles. 

• A quasi-two-body amplitude such as D° — > <pp°, where subsequently 4> — >■ K + K~ 
and p° — > 7r + 7r~. If both intermediate particles are vector mesons than they can be 
in either an S, P, or D wave orbital angular momentum state, which leads to three 
distinct amplitudes. 

• A cascade amplitude such as D° — > i^i(1270) + K~ , where subsequently fCi(1270) + — > 
K*(892)°tc + followed by If (892)° -> K + tc~. If the first and second intermediate 
resonances are both spin one and of opposite parity the second intermediate resonance 
and the pseudoscalar meson can be in be either in an 5* or D wave orbital angular 
momentum state, which leads to two distinct amplitudes. 

In contrast, for a three-body decay only single-resonance amplitudes, with unambiguous 
orbital angular momentum assignment, are possible. 

The total amplitude for the D° — > K + K~tt + tt~ decay, Ado, is modeled as a coherent 
sum over the i intermediate states considered 

A D o{ai,Pj) = y^OjA(Pj) , (1) 

i 

where aj = \a^\e % ^ 1 is a complex factor and Ai(pj) is a parametrization of the intermediate- 
state amplitude. In the likelihood fit the real and imaginary parts of a; are determined, rather 
than | Oj] and 0j. This is because |aj| and 0, are bounded and cyclic variables, respectively, 
which can lead to numerical problems in the fit. The parametrization of a signal- resonance 
amplitude is given by 

AiiPj) = g t (p j )S i (p j )J r l (p j )J r D (p j ), (2) 

where Gi(pj) and S,i{pj) are the lineshape and spin factor for the resonance, respectively. Here 
J^iiPj) and J-'oiPj) are the angular- momentum barrier penetration factors for the resonance 
and D°, respectively. For most resonances the lineshape is parametrized by a relativistic 
Breit-Wigner propagator with a width that depends upon the spin and the daughter mo- 
menta of the resonance j20[. The only exception is the /o(980) — > 7t + tt~ resonance where a 
coupled-channel (Flatte) lineshape [2l| is used. The values of the mass and natural width 
used in the Breit-Wigner propagators are taken from Ref. The parameters used to 

describe the /o(980) resonance are taken from Ref. 0]. Non-resonant states in which there is 
orbital angular momentum among the daughters are modelled as a very broad resonance with 
a mass corresponding to the measured invariant mass and a width that is very much greater 
than the mass of the D° meson; this leads to the spin factor alone altering the distribution 
of the events over phase space. The spin factors are Lorentz-invariant matrix elements that 
describe angular-momentum conservation in the decay and are described in Appendix |A] 
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The functional form of J 7 is that presented by Blatt and Weisskopf in Ref. 22j, which de- 
pends on both the spin and daughter momenta of the resonance. The parametrization of a 
quasi-two-body amplitude or a cascade amplitude is given by 

A = GltoAGfoASifafflfaffiipArDipj), (3) 

where Q\ (pj) and F( (pj) (j = 1, 2) are the lineshape and Blatt- Weisskopf angular-momentum 
barrier penetration factors for the two intermediate resonances that participate in the decay, 
respectively 

In the combined fit to events tagged as either D° and D° it is assumed that there is 
no CP violation in the decay Therefore, the amplitude for D° — > K + K~ix + -k~ , Ajy>(pj), 
is identical to Ano(pj), except that the charges of the daughters are conjugated; this is 
equivalent to the following interchanges of the four-momenta: p\ <H- P2 and p^ 0^4. 

The D° signal PDF Sr>o{cii,Pj), is given by 

<j /_ _ x e(p j )|^ D o(^,p J )| 2 fi 4 fa) U) 
W*>Pi) j e{p )\A D »{a u p 3 )\m A {p 3 )dp ' W 

where e(pj) is the acceptance parametrized in terms of the four-momenta. The D° signal 
PDF, Sjyj, is identical apart from the substitution of Ajys for Ad°- A method that does 
not require explicit evaluation of the functional form of e is used to fit most data sets. This 
method is described in Section IIII B[ 

The background PDF, B(pj), is determined for each data set either from sideband or 
simulated data. The PDF consists of both combinatoric components and those from specific 
resonances, which are added incoherently. The results of this parametrization are given in 
Sec. IIV1 Therefore, the log-likelihood function for a flavor-tagged data set is 

InC = J2 ln [fs{( 1 -^ S Do(a l ,p k j )+uS w (a t ,p l ;)} + (l-fs)B(p k j 

k 

^ In [f s {(l-u)S m (a l ,p k J )+uS D0 (a l ,p k j )} + (1 - fs)B(p k )} , (5) 



where fs is the fractional amount of signal in the data sample, u is the mistag rate, p k 
are the four- momenta of the D° and its daughters for the k th event, and Npo (N-^o) is the 
number of events tagged as D° (D°). 

The signal PDF for CP-even («S + ) and CP-odd (<S_) tagged data is given by: 

s ^ a s = eiPiMApojauPj) ± A w (ai,p j )\ 2 R 4 {pj) 
± *' 1 f e(p j )\A D o(a h p j )±A^(a h p j )\ 2 R i {p j )dp j ' 

Hence, the log-likelihood function fit to CP-tagged data is 

N+ 

ln£ CP = ^ln[/+5 + (a i) p*) + (l-y3-)B + (p*)] 

k 
JV_ 

+ ^ In [fsS_{a h p k ) + (1 - fs)B-(p])] , (7) 
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where N + (AL), fg (fg) and B + (BJ) are the number of tagged decays, fraction of signal, 
and background PDF for the CP-even (CP-odd) eigenstate, respectively. 

The final type of data used in the fit is the CLEO-c 3770 sample of mixed CP tagged 
by K sl 7t + 7t~ . The quantum-correlated amplitude for K + K~tt + tt~ vs. K s 7t + tt~ decays is 
given by: 

[A D «($)A{m 2 _, m 2 + ) - A w (p^)A(ml, ml)] , (8) 

where A(m 2 + , m 2 _) is the amplitude for D° — > K s tt + tc~ as a function of Dalitz plot variables 
m 2 + and ml, which are the invariant-mass squared of the K s n + and K s tt~ pairs. (The 
amplitude for D° — > K s 7i + 7T~ is equal to A{m 2 i : mZ r ) assuming there is no CP violation 



in the decay.) Following Refs. 23|, |24|, dividing the D° — > Kg7i + n^ Dalitz plot into bins 
symmetrically about the line m 2 + = m 2 _ allows the mixed-CP amplitude squared for D° — > 
K + K~7t + tt~ in the m th bin to be written as 



\An(Pj)\ 2 oc \ A D o\ 2 K_ m + \ A w \ 2 K m - 2^/K m K_ m [c m Re(A D oA^o) + s m l^(A D oA^)] , 

(9) 

where K m is the fraction of D° — > Kg7r + ir decays in the m th bin and c m (s m ) is the 
amplitude-weighted average of the cosine (sine) of the strong-phase difference between D° 
and D° decays to Kgn + 7r" within the bin. The pairs of symmetric bins have index m if they 
lie in the region m 2 _ < m 2 + and — m in the region ml > m 2 ^. The values of K m , c m , and s m 
used in this analysis are those measured by the CLEO Collaboration [12'!] . There are several 



binnings presented in Ref . [12j ; this analysis uses the binning in equal intervals of the strong- 



phase difference, which is referred to as the equal A5n binning in Ref. [12] and herein. 



There is an equivalent expression for the amplitude of events tagged by D — > K l h + tx : 

\A' m ( Pj )\ 2 oc \A D o\ 2 K'_ m + \A w \ 2 K' m + 2y/K' m K'_ m [c' m Re(A D oAy + sX^^)] , 

(10) 

where K' m , c' m , and s' m are analogous parameters to those defined for D° — > KgiT + TT decays. 
The value of these parameters used are those for the equal A5d binning reported in Ref. |l2| . 
The signal PDF for Kg7r + 7r~ (fT°7r + 7r~), Sm, is given by: 

sua*,*) = , ! fa)l lE (at ' B) l^ 4fa) . ■ (ii) 



/ e (pj ) | Ahl (a* , pj ) 1 2 P 4 (Pj ) dp 



j 



The log-likelihood function that is maximized for the K L s ir + ir tagged event sample is 



Nrr 



In ^P-mixed = E 111 [fs^^m^ pj) + (1 - /^^CP-ruix^)] 

m=—8,mj^0 (. k 

N' ^ 

+ E ln [fs^'SL^p*) + (1 - f^-^B'cp^UP-)} \ , (12) 

where Nm , f^ p ~ mix ^ anc i B^ p __ mix are the number of tagged decays in each bin, fraction of 
signal, and background PDF for the Kgir + 7T~ (AT°7r + 7r~) tagged events, respectively. 

The combined log-likelihood function to be maximized is the sum of the expression in 
Eq. §5§ for each flavor-tagged data set, plus the log-likelihood functions given in Eqs. ([7]) 
and (j!2p . While performing the fit, one intermediate resonance component has dj fixed to 
unity such that the amplitude and phase of the other components are determined relative 
to it. 
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B. Efficiency parametrization and normalization 



The efficiency parametrization and normalization methods are presented together in this 
section as they are closely linked. Consider first the logarithm of the flavor-tagged D° signal 
PDF given by Eq. (0}: 



In S D o (cii , pj ) = In e (pj ) + In | A D o (a* , Pj ) | 2 + In R A (pj ) - In 



e(pj) \A D °{aiiPj)\ 2 R±{Pj)dpj 



(13) 

Only the second and last terms depend on the amplitude model parameters, which means 
that only these need to be computed while performing the fit. Therefore, the acceptance 
function, e{pj), is only required to compute the normalization integral. The normalization 
integral is evaluated by a Monte Carlo integration method, as described below, which does 
not require an analytic form of e(pj). Not requiring an analytic form is desirable, given such 
acceptance functions can be difficult to parametrize even over a three-body phase space, due 
to the rapid change in efficiency at the edge of phase space. These problems are compounded 
for four-body decays due to the higher dimensionality of the phase space. 

However, the treatment is less straightforward in the presence of background. The com- 
bined PDF, ignoring the mistag rate, is 

f s S I (1 f s )B = fs e (Pi^ AD °( ai >P^ 2Ri (P^ [_(! f s \ gfeWgij ( 14 ) 
S Je(p j ) \A D o(ai, pj^R^dpj s J B{p j )R^{p j )dp j ' 

where B(s) is a function that describes the background distribution relative to phase space. 
Given that e no longer factorizes from the second term it has to be computed for each event 
to minimize the log-likelihood function. The combined PDF becomes 



f s S + (l-fs)B = e(p j )R 4 {p j 



\A D o(ai,pj)\ : 

Js- 



' fe(pj)\ A D o (a, , pj ) | 2 Ra (pj ) dp 



j 



+(1 - fs 



f4Pj) B e(Pj)R4(Pj)dPj 



(15) 



where B e (pj) is defined to be the background distribution relative to the acceptance corrected 
phase space [^(Pj)R^(Pj)]- Now that e factorizes from the complete PDF, it only has to be 
computed as part of the normalization. 

The normalization integrals are determined by a Monte Carlo technique. The simulated 
events are generated according to the distribution 

\A D0 (ar, Pj )\ 2 R 4 ( Pj ) , 

where af en are a fixed set of parameters. Using this set of simulated events the normalization 
integral can be approximated by 



/ eip^An^p^R^) ~ J- g e(p)) I , (16) 

J iv gen k |^D°l u i iPj}\ 



where iV gen is the number of simulated events. 



Furthermore, the effect of the acceptance function can be incorporated in determining this 
integral by summing over events, that once passed through the CLEO detector simulation, 
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satisfy the selection criteria described in Sec. [IT] for a given data set. This is equivalent to 
generating events according to the PDF: 

Therefore, the normalization integral is given by 



/ 



< MA < ^ 1 V* \^no(a^)\ 2 

efe)|^oKp,)l RAPi) ^L ^^^p . ( 17 ) 



where iV se i is the number of simulated events selected. The background samples are also fit 
using these normalization events so that the function B e is determined. 

The number of integration events is chosen to ensure the uncertainty on the integral 
is less than 0.3%. Two separate samples of one million selected events each are used to 
perform the integration of the fits to CLEO III and CLEO-c data. It is no longer possible to 
generate an additional signal simulation sample of this size for the CLEO II. V data. Given 
the similarity between the CLEO II. V and CLEO III detectors the CLEO III integration 
events are used. To account for the small differences in particle identification and tracking 
performance between CLEO II. V and CLEO III, the integration events are reweighted by 

n £ H.v(|Pjl) / 18 n 

where en.v(IPjl) [ e ni(|Pil)] is the efficiency as a function of the momentum of the daughters 
in the laboratory frame, |pj|, for CLEO II. V [CLEO III]. The ratios en.v/ e m are computed 
from simulation. The mean value of the weight applied is 0.87 with an R.M.S. of 0.04. Such 
a reweighting does not account for correlations among the daughters and the dependence 
of the acceptance on other variables; therefore, a conservative systematic uncertainty is 
assigned to account for this approximation, which is discussed further in Sec. IIVDI 



C. Goodness-of-fit 

In order to quantify the quality of a given fit a x 2 value is computed. The four-body 
phase space can be described completely by any five invariant-mass-squared variables s^- = 
(Pi + Pj) 2 an d Sijk = (pi + Pj +Pk) 2 - Therefore, the events are binned in terms of Si 2 , s 12 3, 
S23, S234, and S34 to compute the x 2 ■ Initially, the phase space is divided into equal bins. 
At least fifty events are required in each bin so that the x 2 calculation is robust. Therefore, 
after the initial equal division of the phase space, bins are merged until they satisfy the 
minimum number of events criterion. 



The x 2 is given by 



E 



where N p and Np Xp (di) are the observed and expected number of events per bin, respectively, 
and n is the number of bins. In the general case the value of N!: xp (ai) is given by 

N^{oi) =N f [fsS D0 ( ai ) + (1 - f s )B] d Pj , (20) 

J bin p 
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where N is the total number of events in a particular sample and the integral is over the p 1 
bin. The Monte Carlo integration events are used to compute this integral such that 



iVf P (a 4 ) 



iVsel 



N 




fslApoja^ 2 + (1 - f s )B e 
\A D o{af en )\ 2 



(21) 



where N^ el is the number of Monte Carlo integration events in the p th bin. 
The number of degrees of freedom, u, is given by 



v = 



(n — 1) — n 



par 



(22) 



where n par is the number of free parameters in the amplitude or background model being fit. 
The number of bins is reduced by one because the number of expected events in the final 
bin considered is determined by the overall normalization and as such does not represent a 
degree of freedom. The compatibility of the combined fit hypothesis to each individual data 
set is estimated by determining the x 2 P er bin for that data set. The CP-tagged samples 
are not used in computing the \ 2 because of the limited statistics in the four sub-samples. 

IV. RESULTS 

The results of the amplitude analysis of D° — > K + K~7t + tt~ are presented in this section. 
The fits to the sideband samples used to parameterize the background in each data set 
are presented in Sec. IIV Al The method used to arrive at the components in the baseline 
model and the results are presented in Sec. IIVBI Robustness tests of the fitting method and 
cross checks of the final result are given in Sec. IIV CI The procedures to evaluate systematic 
uncertainties are outlined in Sec. IIVDI A search for CP violation in the decay is presented 



A. Background parameterization 

The fraction of signal in each data set considered in the amplitude fit has been estimated 
for the various selections and is given in Table [TTJ As well as the different relative amount of 
background in each data set, the composition of the background is also different due to the 
variation in the value of y/s, the detector configuration and tagging method. In addition, 
the type of backgrounds must be classified into those that form candidates that are peaking 
or non-peaking at the nominal D° meson mass. 

Non-peaking backgrounds are random combinations of four particles that do not originate 
from the same D° decay, but some may form a resonance such as a or p. Separate sideband 
samples have been selected as described in Sec. [Til which are fit to determine the non-peaking 
background PDFs for each data set. The model used to describe the background is an 
incoherent sum of resonances along with a non-resonant component, which is a constant, 
such that 



where b nr and 6j are real parameters determined by the fit. Various combinations of am- 
plitudes are tested. The combination with the lowest x*/v, which does not contain any 



in Sec. HVEl 




(23) 



17 



components that contribute less than 0.5% to the total, is selected. Table [TTT] gives the 
fractional contributions of the different components for each flavor-tagged data set. All data 
sets have a significant non-resonant component in the background but the resonances that 
contribute vary significantly among the different data sets. 

TABLE III: Fractional contribution of each component to the background for each flavor-tagged 
data set. The x 2 /v of the background fit is also given. 







CLEO II. V 


CLEO III 


CLEO-c 3770 


CLEO-c 4170 


#1(1270)- 


(#J(1430)°7T- 


)#+ 


— 


0.019±0.064 


— 


#i(1270)+ 


(#*°7T+)#- 


0.053±0.015 








#1(1270)" 


(#*°7T-)#+ 


0.006±0.015 


0.032±0.015 


0.033± 0.021 




#i(1270)+ 


(p#+)#- 






0.005±0.013 




#i(1400)+ 


-(#*°7T+)#- 




0.055±0.014 


0.015±0.023 










0.079±0.011 


0.143±0.022 


0.102±0.013 


#*°#*° 




0.007±0.012 


0.045±0.014 


0.101±0.019 


0.010±0.017 


#*°#+7T- 










0.032±0.023 


/o(980)#+ 


#- 


0.033±0.047 


0.128±0.049 


0.293±0.118 




pK+K- 










0.243±0.034 


K*°K-tt+ 








0.017±0.035 


0.098±0.025 


Non-resonant 


0.899±0.048 


0.661±0.046 


0.373±0.103 


0.516±0.078 


X 2 /v 




1.32 


1.20 


1.17 


2.13 



The only significant peaking background comes from D° — > Kg(7T + 7i^)K + K^ decays in 
the CLEO-c datasets. The larger average momentum of the K° s mesons in the CLEO II. V 
and CLEO III data sets leads to a significant displacement of most Kg decay vertices from 
the interaction point, which allows the efficient rejection of this background. The lower 
average momenta at CLEO-c means that the such a separation is less effective. The fraction 
of this component is left as a free parameter in the fit to data. A four-body model of 
the distribution of the Kg(7T + 7T~)K + K~ events over the kinematic variables is obtained by 
fitting events that fail the K$ veto. The resonant components of the model considered are 
a subset of those reported in Ref. [25[ and these are fit coherently to this sample. Only 
those with a significant non-zero contribution are retained. The model parameters found 
are then used in the fits to the signal sample data to yield the background fractions reported 
in Sec. [ITBT1 and IITB21 

The non-peaking background distribution for the CLEO-c CP-tagged data is estimated 
from a fit to a combination of generic simulation events and data sidebands, because of the 
limited statistics. There is also a peaking D° — > Kg(7r + 7r~)K + K~ contribution, which is 
estimated from the simulation alone. Due to the reliance of these background estimates 
on simulated events the strategy to determine the systematic uncertainty related to the 
background in the CP-tagged sample is more conservative than the flavor-tagged samples; 
this is discussed further in Sec. HVDl 
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B. Model selection and results 



There are many possible amplitudes which can contribute to the D° — > K + K^-k + -k + 
decay, therefore a strategy to determine the best combination is defined. From inspection 
of the invariant- mass squared projections (see Figs. and ED there is clear evidence for 
intermediate <fi — > K + K~, K*° — > K~tt + , p° — > 7r + 7r~ and a broad high-mass kaon decaying 
to K-jv+Tv- (^(1270)", K^UOO)-, K*(U10)-, K*(U30)-, or K*(1680)~). Therefore, all 
models considered contain at least one amplitude with one of these intermediate resonances. 
In addition, for states that are not self conjugate, the conjugate amplitude is always included 
as well; for example, if Ki(1270)~ K + is a component in the model so is Ki(l270) + K~ . 

Models containing seven components are tested and fit to the combined data set. The 
five models with the smallest x 2 l v are considered further. An additional component is 
then added, which either contains one of the principal resonances discussed above or is non- 
resonant. All models, including those from the previous iteration, are compared and the 
best five are retained. This process continues until the five models with the lowest x 2 1 v are 
the same as those in the previous iteration. 

At this point amplitudes containing intermediate resonances not already considered, such 
as u), /o(980), and /2(1270), are added. The best fifty models are then retained. Any models 
that contain components contributing less than 5% are then simplified by removing these 
components and the revised model is fit to the data. Testing these simplified models ensures 
that the improvement in the y 2 / v by including these small components is significant. Finally, 
all models that have been tested are ranked according to their x 2 1 v - The different amplitudes 
that have been included in the models tested are listed in Appendix [B] 

The components of the model with the lowest x 2 l v are shown in Table IIVI For some 
amplitudes, pairs of particles do not decay via a resonance but are in a state of relative 
orbital angular momentum (L); each such pair is surrounded by curly brackets with an S or 
P subscript, indicating an L = or L = 1 state, respectively. The relative orbital angular 
momentum state - S, P, or D wave - of intermediate resonances and pairs of particles is also 
given for an amplitude if more than one is possible. There are other models with similar 
X 2 l v ] the principal variations are different or additional angular momentum states for the 
07r + 7r~ decay. More information on these alternative models can be found in Appendix 
The x 2 l v f° r the combined fit and the x 2 /bin for each flavor-tagged data set are presented 
in Tabled 

The best fit projected on to the distributions of Sij and Sijk for the combined flavor-tagged 
data set are shown in Figs. |5] and [61 repsectively. Reasonable agreement is seen between the 
data and the fit for most distributions. An exception is the narrow peak in the background 
S34 distribution, which is due to the peaking K®K + K~ background found in the CLEO-c 
data. The best fit underestimates the data in this region; therefore, an additional systematic 
uncertainty, described in Sec. UVD[ is assigned to account for this discrepancy. 

The best fit real and imaginary components of a, are given in Table IIVI The statistical 



correlations among the real and imaginary components of <2j are given in Ref. [26(. The 
values of the magnitude and phase of the amplitude derived from the fitted parameters are 
given in Table I VI I along with the associated statistical and systematic uncertainties. The fit 
fraction, which is defined as 

/ \M 2 dpj ( 
f\A D o\*d Pj > [ ^ } 

is also given for each component. The fit fraction indicates the relative contribution of 
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TABLE IV: Real and imaginary parts of for combined fit to all data. Only the statistical 
uncertainties are given. The daughters of the -RTi(1270) =l: are assumed to be in an S-wave state. 



Amplitude 


Re(aj) 


Im(aj) 


i^i(1270) + (K*°vr + )K- 


1.0 


0.0 


i^i(1270)-(K*°7r-)K+ 


0.16 ±0.08 


-0.31 ±0.06 


Ki(l27Q)+lp K+)K- 


4.07 ±0.64 


4.22 ±0.89 


Ki(1270)~ (p°K-)K + 


6.90 ± 0.59 


0.20 ± 1.10 


K*(U10) + (K*°-k + )K- 


4.62 ±0.56 


-4.10 ±0.72 


K*(U10)-(K*°ir-)K + 


2.61 ±0.93 


-6.25 ±0.59 


K*°K*° S wave 


0.32 ± 0.04 


0.13 ±0.04 


4>p° S wave 


-0.32 ±0.15 


0.99 ±0.09 


4>p D wave 


0.20 ±0.33 


-1.43 ±0.18 




-1.70 ±0.83 


-5.93 ±0.48 


{K-n+} p {K+n-} s 


82.6 ±6.8 


11.4 ±9.7 



TABLE V: x 2 l v f° r the combined fit or x 2 / n f° r each flavor-tagged data set. For the combined fit 
there are 22 free parameters. 



Data set 


X 2 /u or x 2 /n 


v or 


Combined 


1.63 


113 


CLEO-c 3770 


1.29 


55 


CLEO-c 4170 


1.20 


26 


CLEO III 


1.54 


49 


CLEO II. V 


2.00 


6 



a component to the total branching fraction. The individual component fractions do not 
have to sum to 100% due to interference. (Information about the interference among the 



amplitudes is given in Ref. |26f] . ) For the best fit the sum of the fit fractions is (96.7 ± 
2.6 ± 9.8)% where the first uncertainty is statistical and the second is systematic. The 
dominant intermediate state is the (ftp quasi-two-body decay. There are three other decay 
modes which contribute over 10% to the total branching fraction: Ki(1270) + K~, 07r + 7r~, 
and the non-resonant (K~7r + )p(K + 7r~)s- 

Detailed comparison of these results to those presented previously [H, @] is not straight- 
forward given the lack of flavor tagging in both analyses and the absence of spin factors in 
the E791 study. However, there is agreement with the previous findings in the presence of 
a significant 0p° contribution. Further, a significant 07r + 7r~ is observed by the E791 collab- 
oration [l[. Significant contributions from both the <pp and 07r + 7r~ modes are anticipated 
because there is only one singly-Cabibbo-suppressed diagram that contributes to the rate 
[lj]. However, for other modes such as K*°K + n~ there are two leading-order amplitudes of 
opposite sign, which result in a suppression of the decay rate. In addition, there is evidence 
of a large Ki(\270) + K~ contribution, which was also observed by the FOCUS collaboration 
(2)]. However, the expectation that K\ (1270) mesons with the same charge as the W boson in 
the decay will be the dominant amplitude |2| is not observed. This may be a consequence of 
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TABLE VI: Modulus, phase, and fit fraction for each component of the baseline model. The first 
uncertainty is statistical and the second is systematic. 



Component 








\ai\ 


#1 (1270)+ (ir^V 


^)K- 






1.0 


Ki(1270) - (K*°ir- 


-)K+ 


0.35 


± 


0.06 ± 0.03 


#1(1270)+ (j>° K + )K- 


5.86 


± 


0.77 ±2.03 


K l (1270)-{p°K- 


)K+ 


6.90 


± 


0.59 ± 3.07 


K*(U10) + (K*°it- 


+)K~ 


6.18 


± 


0.64 ±0.75 


K*(UW)-(K* n 


-)K+ 


6.78 


± 


0.65 ± 1.25 


K*°K*° S wave 




0.34 


± 


0.04 ±0.14 


4>p° S wave 




1.04 


± 


0.10 ±0.31 


4>p° D wave 




1.44 


± 


0.19 ±0.38 


0{7T+7r-} 5 




6.17 


± 


0.52 ± 1.58 


{K-TT + } p {K + TT- 


■h 


83.4 


± 


6.8 ±29.3 



4>i (rad) Fit Fraction (%) 



0.0 


7.3 ±0. 


.8 ±1.9 


1.10 ±0.22 ±0.23 


0.9 ±0. 


,3 ±0.4 


0.80 ±0.13 ±0.08 


4.7 ±0. 


,7 ±0.8 


0.03 ±0.16 ±0.23 


6.0 ±0. 


,8 ±0.6 


0.73 ±0.11 ±0.33 


4.2 ±0. 


,7 ±0.8 


1.18 ±0.13 ±0.48 


4.7 ±0. 


.7 ±0.7 


0.39 ±0.12 ±0.18 


6.1 ±0. 


,8 ±0.9 


1.89 ±0.14 ±0.35 


38.3 ± 2. 


,5 ±3.8 


1.43 ±0.22 ±0.48 


3.4 ±0. 


,7 ±0.6 


1.85 ±0.13 ±0.37 


10.3 ±1. 


,0±0.8 


0.14 ±0.12 ±0.28 


10.9 ±1. 


.2 ±1.7 



final-state interactions playing a significant role. Further evidence of final-state interactions 
is given by the presence of a statistically significant K*°K*° contribution, because in the 
SU(3)-flavor limit, the two ^-exchange amplitudes contributing to this final state cancel 

Q. 

There are two significant differences from the FOCUS model. Firstly, no significant 
/o(980)7r+7r~ contribution is found. Secondly, a non-resonant contribution with angular 
momentum structure is required to fit the data. (If there is no angular momentum among 
the particles in the non-resonant component the x 2 ' l v increases to 2.1.) 

It is of interest to assess the impact of the CP-tagged data given that the quantum- 
correlated states should provide additional information about 0j compared to the flavor- 
tagged data. Given the limited number of events in the sample it is not possible to fit the 
CP-tagged data alone. However, the CLEO II. V flavor-tagged data sample has a similar 
size and purity to the CLEO-c CP-tagged data sample. Therefore, the statistical impact 
of these two samples is compared by determining how much the mean relative statistical 
uncertainty on Re(aj) and Im(aj) changes when either the CP-tagged CLEO-c data set or 
the flavor-tagged CLEO II. V data set is excluded from the combined fit. The mean relative 
statistical uncertainty of the increases by 12.3% when the CP-tagged data set is removed 
compared to 7.5% when the CLEO II. V data set is removed, indicating that a small sample 
CP-tagged data is more powerful than an additional sample of flavor-tagged events of similar 
size. 



21 



4280112-013 



> 

CD 

CD 

CO 

o 

CD 

t/> 

CD 
> 
LU 



> 
CD 
O 

CM 
O 
CD 

CO 

H — ' 

CD 
> 
LU 



600 
500 
400 
300 
200 
100 




350 
300 
250 
200 
150 
100 
50 




: i 




(a) 


1 1.2 1.4 1.6 1.8 2 


2.2 2.4 2.6 




S 12 (GeV 2 /c 4 ) 








(c) 




jj \ 












r \ 











0.4 0.6 



0.8 1 1.2 

S 14 (GeV 2 /c 4 ) 



> 

CD 
O 

CM 
O 
CD 

03 

H — I 

c 

CD 
> 
LU 



1.4 1.6 



220 






(e) 


200 






180 








160 








140 








120 








100 


[ / ^ 






80 








60 








40 








20 











0.4 0.6 0.8 


1 


1.2 1.4 1.6 



> 

CD 

h- 

CM 
O 

CD 



c 

CD 
> 
LU 



> 
CD 
O 

r- 

CM 
O 
CD 

i2 
c 

CD 
> 
LU 



> 

CD 

o 

o 

CD 



CD 
> 
LU 



250 
200 
150 
100 
50 




300 





(b) 


\ I \ 




L / T*t 








0.4 0.6 0.8 1 1.2 


1.4 1.6 


S 13 (GeV 2 /c 4 ) 






S 23 (GeV 2 /c 4 ) 




S 24 (GeV 2 /c 4 ) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

S 34 (GeV 2 /c 4 ) 



FIG. 5: The (a) S12, (b) S13, (c) su, (d) 523, (e) S24, and (f) S34 projections for all flavor-tagged 
data (points with error bars) with the best fit (solid line) superimposed. The indices correspond 
to K + = 1, K~ = 2, 7r + = 3, and ir~ = 4. The contributions from mistag (filled region) and 
background plus mistag (dashed line) are also shown. 



22 



4280112-014 




S 134 (GeV 2 /c«) S 234 (GeV 2 /c 4 ) 



FIG. 6: Distribution of for the three-body invariant-mass observables: The (a) S123, (b) S124, (c) 
S134, and (d) S234 projections for all flavor-tagged data (points with error bars) with the best 
fit (solid line) superimposed. The contributions from mistag (filled region) and background plus 
mistag (dashed line) are also shown. 
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TABLE VII: Fit fractions for the CLEO II.V/CLEO III and CLEO-c 3770/4170 data sets fit 
separately and the difference between the two results, normalized by the uncorrelated uncertainty. 



Component 




Fit fract 


ion (%) 


Difference (cr) 






CLEO-II.V/III 


CLEO-c 




Ki(1270) + (ir"°7r 


-)K- 


4.4 ±0.9 ±0.8 


9.8 ± 1.2 ±4.6 


1.1 


Ki(1270)~ (K*°ir- 


-)K+ 


3.6 ± 0.9 ± 1.1 


0.2 ±0.2 ±2.3 


1.2 


Ki(1270) + (p°K + )K- 


2.6 ±0.8 ±0.5 


8.5 ±1.4 ±4.2 


1.3 


K l (1270)-{p°K- 


)K+ 


7.9 ±1.2 ±0.8 


3.5 ± 1.1 ± 3.4 


1.2 


(1410)+ (if 


^)K- 


4.5 ±0.9 ±0.6 


4.5 ± 1.1 ± 1.2 


0.0 


K*(UW)-(K*°tt- 


-)K+ 


5.5 ± 1.0 ±0.8 


4.6 ±0.9 ±1.1 


0.4 


K*°K*° S wave 




7.5 ± 1.8 ± 1.8 


5.2 ± 1.0 ± 2.0 


0.7 


<f>p° S wave 




39.8 ±2.7 ±1.4 


36.9 ±3.2 ±4.0 


0.5 


<pp° D wave 




4.7 ± 1.0 ± 1.1 


1.7 ±0.8 ±2.3 


1.0 


<t>{ir + -K-} s 




8.3 ± 1.2 ±0.5 


10.9 ±1.7 ±3.0 


0.7 


{K--k + } p {K + tt- 




14.7 ± 1.8 ±4.3 


7.8 ± 1.4 ± 5.3 


1.1 



C. Robustness tests 

Tests of the result are made by dividing the data into different subsamples. Firstly, 
separate fits to the data collected at center-of-mass energies around the T resonances (CLEO 
II. V and CLEO III) and data collected at CLEO-c are made. This tests the compatibility 
of the results using data produced at different center-of-mass energies and selected with 
different tagging techniques. The resulting fit fractions are compared to one another in 
Table IVHI The number of standard deviation difference between the two results takes into 
account the uncorrelated systematic uncertainties between the two data samples. There is 
good agreement between the two samples. 

The fitter is also tested on an ensemble of simulated data sets to identify any bias in 
the fit and determine the reliability of the statistical uncertainties returned by the fit. The 
ensemble is 200 sets of 4000 events generated with the CLEO-c simulation, which are fit 
individually. The distribution of aj for the ensemble of experiments is shown in Fig. 
the values are seen to be scattered about the generated values with no significant biases. 
The distribution of the pull, defined as the difference between the fitted and generated 
parameters divided by the uncertainty on the parameter returned by the fit, is formed for 
the real and imaginary part of a, for the ensemble of simulation experiments. If the fit is 
unbiased and the uncertainties are correctly determined the pull distribution will be normal. 
This is tested by fitting each pull distribution with a Gaussian function. The width and 
mean found by the fit to the pull distributions are given in Table [VTTT1 The ensemble study is 
performed with the fit using the generated and reconstructed four-momenta of the particles. 
As the resolution is neglected in the fit the difference between the results obtained with 
the generated and reconstructed four-momenta allows the systematic uncertainty related to 
resolution to be determined. Small biases are observed but only one is greater than three 
standard deviations from zero. The widths of the pull distributions are all compatible with 
unity indicating the uncertainties are correctly evaluated. The small bias is accounted for in 
evaluating the systematic uncertainties, which is described in Sec. IIVDI The average \ 2 l v 
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FIG. 7: Fitted <Zj (points) distribution for the ensemble test of 200 simulated data 
sets. Also shown is the generated value of (asterisk). The labels correspond 
to: (a) jjfi(1270)- {K*^tt~)K + , (b) ifi(1270) + (p° K + )K~ , (c) K^mO)- (p K~)K + (d) 
K*(UW) + (K*°tt + )K-, (e) K*(UW)-(K*®7r-)K + , (f) K*°K*° S wave, (g) <j>p Q S wave, (h) cpp 
D wave, (i) (p{ir + TT^} s , and (j) {K~ir + } p {K + n~} s . The values of Re(dj) and Im(aj) are scaled 
by 0.1 for the {K~ir + } p {K + ir~} s amplitude. 



for the ensemble of fits is 0.96 ±0.01, which is a further indication that the fitting algorithm 
is well behaved. 



D. Systematic uncertainties 

The systematic uncertainties on the results fall into several categories: amplitude model 
assumptions, parameterization of the background, modeling of experimental effects, and the 
fitter performance. Each of these categories is discussed below. The resulting systematic 
uncertainties for the amplitudes, phases, and fit fractions due to each contribution are given 
in Tables ED and [3 

Three assumptions of the amplitude model are tested: the mass and width of resonances, 
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TABLE VIII: Mean (p) and width (a) of the pull distributions from simulated data for all of the 
fitted parameters using either the generated or reconstructed four-momenta. 



Parameter Generated pi Reconstructed pi 







A* 








a 














a 




Ki(1270)~(K* tt-)K+ Re(ai) 


-0 


43 ±0 


13 


1 


02 


±0 


14 


-0 


36 ±0 


12 





99 


±0 


13 


Ki(1270)-(K \-)K+ haifli) 





03 ±0 


13 


1 


14 


±0 


12 


-0 


10 ±0 


14 


1 


13 


±0 


13 


K"i(1270)+ (p° 'K + )K~ Re(oj) 





08 ± 


12 


1 


06 


± 


12 





16 ± 


19 


1 


35 


± 


19 


i / r\f~7f~\\ -4- / TV-4- \ TV* — T f \ 

Ai(1270) + (/9 u A + )A Im(aj) 


— 


15 ± 


09 





85 


± 


07 


— 


29 ± 


10 





83 


± 


10 


Ai(1270) (p K )K + Re(oj) 


— 


01 ± 


11 





88 


± 


10 


— 


03 ± 


12 


1 


01 


± 


11 


Ki (1270) (p^K )K + Im(aj) 


-0 


24 ± 


14 


1 


17 


± 


14 


-0 


22 ± 


12 


1 


08 


± 


12 


/t (1410) (iv 7T ), iv Re(aj) 


— 


32 ± 


13 


1 


11 


± 


12 


— 


41 ± 


11 





98 


± 


09 


K (1410 y(K u Tr + ) K Im(aj) 





08 ± 


14 


1 


08 


± 


12 





11 ± 


11 





89 


± 


10 


K*(1410) _ (i^ tt-)K + Re(aj) 


-0 


05 ±0 


11 





91 


±0 


10 


-0 


08 ± 


13 





96 


±0 


11 


— *n 

A"* (1410)" (A 7T — ), .ff Im(ai) 


-0 


35 ±0 


11 





95 


±0 


12 





23 ±0 


13 


1 


06 


±0 


12 


K*°K S wave Re(oi) 





28 ±0 


11 





94 


±0 


09 





21 ±0 


14 


1 


13 


±0 


12 


K"* ^* S wave Im(oi) 


-0 


15 ±0 


11 





99 


±0 


10 


-0 


02 ±0 


10 





83 


±0 


09 


S wave Re(aj) 


-0 


29 ±0 


11 





87 


±0 


15 


-0 


26 ±0 


10 





90 


±0 


12 


<fip° S wave Im(cij) 


-0 


02 ±0 


10 





85 


±0 


09 





17 ±0 


10 





86 


±0 


10 


D wave Re(oj) 





20 ±0 


13 


1 


05 


±0 


11 





31 ±0 


13 


1 


11 


±0 


12 


(ftp D wave Im(oj) 


-0 


06 ±0 


11 





96 


±0 


09 





01 ±0 


12 





91 


±0 


10 


(f){ir + ir~} s Re(oj) 


-0 


20 ±0 


14 


1 


14 


±0 


12 


-0 


13 ±0 


14 


1 


16 


±0 


12 


0{7r + 7r-} 5 Im(a;) 


-0 


16 ±0 


09 





85 


±0 


10 


-0 


24 ±0 


13 


1 


08 


±0 


12 


{A"-7r+} P {A"+vr-} s Re(ai) 





17 ±0 


10 





87 


±0 


10 





01 ±0 


12 





98 


±0 


11 


{A-vr+} p {A-+vr-} s Im(a 4 ) 





03 ±0 


16 


1 


15 


±0 


21 


-0 


07 ±0 


12 





86 


±0 


11 



the barrier-penetration factors, and the absence of quantum correlations in the modeling of 
CLEO-c flavor tagged data. The mass and width assumed for the resonances in the model, 
Ai(1270) + , K*°, A"(1410) + , <fi, and p°, are taken from Ref. Q; these are varied by the quoted 
uncertainties to determine the related shift of the fit parameters. These shifts are added in 
quadrature to obtain the total uncertainty related to the mass and width parameters. The 
uncertainty from this source is of the same size or larger than the statistical error for most 
parameters. The uncertainty on the A'i(1270) + mass and width dominates. Determining 
the A'i(1270) + mass and width from data does not improve the overall uncertainty. 

The amplitude model includes spin- and momentum-dependent Blatt-Weisskopf orbital- 
angular-momentum barrier penetration factors. These factors are set to unity and the fit 
repeated. The resulting shift in the fitted parameter values and fit fractions is taken as the 
systematic uncertainty. 

Quantum correlations in the CLEO-c flavor-tagged data are ignored in the fit. Since 
many different final states containing the tagging kaon are summed over in the analysis 
the effect of correlations is diluted. However, an alternative signal PDF is tested, which 
accounts for the correlations at the cost of two additional parameters (see Appendix iDl). 
The resulting changes in the central values of the fitted parameters and fit fractions, taken 
as the systematic uncertainties from this source, are found to be less than half a statistical 
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TABLE IX: Systematic uncertainties on \ai\ and (pi in units of statistical standard deviations (a). 
The different contributions are: (I) mass and width of resonances; (II) Blatt-Weisskopf penetra- 
tion factors; (III) quantum correlations; (IV) background fractions; (V) flavor-tagged background 
parameterization; (VI) CP-tagged background parameterization ; (VII) K$K + K~ background; 
(VIII) acceptance; (IX) resolution; (X) mistag rate; and (XI) fitter bias. 



Parameter Source (a) Total (a) 









I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 






i4Ti(1270)-(if*°7r- 


-)K+ \ ai \ 





03 





17 





01 





31 





15 





08 





01 





23 





00 





03 





24 





52 


Ki(1270)-(i^* 7r- 


-)K+ & 





83 





30 





04 





0G 





18 





02 





02 





30 





00 





04 





40 


1 


05 


K l {mQ) + { P °K+ 


)K~ H 


2 


52 





23 





17 





0G 





30 





03 





04 





43 





31 





01 





07 


2 


Gl 


K 1 {mQ)+[p Q K+ 


)K~ & 





46 





14 





11 





13 





17 





0G 





06 





07 





14 





02 





1G 





58 


i^i(1270)-(/J°iC- 


)K+ \ai\ 


5 


10 





30 





04 





22 





21 





01 





02 





76 





05 





02 





01 


5 


17 


K 1 (127Q)-(p°K- 


)K+ fa 


1 


27 





42 





07 





11 





07 





04 





06 





51 





22 





04 





13 


1 


47 


K*(U10)+(K*°tt 




1 


04 





31 





18 





06 





11 





06 





06 





27 





18 





02 





27 


1 


18 


K*(U10) + (K* o t: 


f )K~ & 


3 


04 





47 





09 





11 





22 





09 





04 





29 





15 





04 





11 


3 


10 


K*(U10)-(K*°n 


-)K+ \oi 


1 


66 





7G 





07 





07 





22 





01 





02 





45 





15 





02 





2G 


1 


92 


K*(U10)-(K*°n 


-)K+ & 


3 


4G 


1 


04 





13 





23 





22 





11 





01 





23 





12 





01 





13 


3 


G4 


K*°K*° S wave | 




3 


86 





27 





23 





20 





17 





07 





07 





55 





00 





01 





19 


3 


93 


K*°K*° S wave </ 


H 


1 


20 





G4 





01 





22 





2G 





11 





02 





37 





00 





02 





21 


1 


47 


(f>p° S wave |dj| 




3 


04 





1G 





22 





23 





1G 





07 





06 





43 





00 





01 





1G 


3 


10 


(f>p° S wave 4>i 




2 


42 





47 





24 





14 





13 





15 





12 





35 





00 





02 





33 


2 


54 


(f>p° D wave \cii\ 




1 


89 





13 





27 





19 





14 





03 





08 





59 





07 





01 





11 


2 


02 


<fip° D wave 4>i 




2 


05 





35 





15 





10 





11 





12 





06 





27 





24 





02 





19 


2 


14 


(f){TT+TT-} s \<Xi\ 




3 


01 





02 





07 





18 





20 





02 





09 





37 





11 





01 





28 


3 


0G 


4>{ir+ir-} s fa 




2 


64 





G8 





19 





09 





12 





16 





12 





38 





02 





03 





1G 


2 


77 


{K-TT + }p{K + TT- 


"Is M 


4 


25 





13 





02 





10 





30 





06 





05 





43 





00 





02 





20 


4 


29 


{K-tt+} p {K+n- 


"Is & 


2 


14 


1 


04 





08 





06 





1G 





11 





06 





17 





00 





02 





05 


2 


39 



standard deviation in all cases. 

Both the level and shape of the background are considered when evaluating the systematic 
uncertainty. The fraction of non-peaking background in each flavor-tagged sample is esti- 
mated from data. The value of fs is varied by its statistical uncertainty for each data sample 
in turn and the shifts in the results are added in quadrature to estimate the systematic error 
from this source. The components in the background model for the flavor-tagged data are 
changed such that those that are not statistically significant, defined as those contribut- 
ing a fraction less than three standard deviations from zero, are removed from the model 
and the fit repeated. Each such component is removed in turn and the individual shifts 
are summed in quadrature to obtain the total uncertainty. In the case of the CP-tagged 
data where the background model is derived from data and simulation the fit is repeated 
ignoring the background. The shifts in the fit parameters with respect to the nominal fit 
are conservatively taken as the systematic uncertainty. In Sec. IIVBI it is noted that the fit 
underestimates the level of the K$K + K~ background. To assess the systematic uncertainty 
related to this discrepancy the fit is repeated with the fraction of the KgK + K~ background 
fixed to double that found in the data, which leads to reasonable agreement with data in the 
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TABLE X: Systematic uncertainties on the fit fraction in units of statistical standard deviations (<r). 
The different contributions are: (I) mass and width of resonances; (II) Blatt-Weisskopf penetration 
factors; (III) quantum correlations; (IV) background fractions; (V) flavor-tagged background pa- 
rameterization; (VI) CP-tagged background parameterization; (VI) KgK + K~ background; (VIII) 
acceptance; (IX) resolution; (X) mistag rate; and (XI) fitter bias. 

Fit fraction Source (a) Total (a) 

I II III IV V VI VII VIII IX X XI 

K 1 (1270)+(K* o n + )K- 2.23 0.57 0.16 0.25 0.27 0.06 0.05 0.57 0.00 0.01 0.05 2A0 

ifi(1270)~ {K*\-)K + 1.13 0.41 0.05 0.28 0.11 0.09 0.01 0.15 0.23 0.04 0.28 1.30 

K 1 (1270)+(p a K+)K- 0.81 0.30 0.13 0.09 0.32 0.01 0.03 0.34 0.08 0.02 0.03 1.00 

K 1 (1270)-(p°K-)K+ 0.41 0.32 0.08 0.13 0.09 0.06 0.04 0.51 0.09 0.03 0.05 0.75 

K*(U10)+(K*°Tr + )K~ 0.84 0.16 0.14 0.14 0.28 0.05 0.05 0.33 0.26 0.02 0.35 1.08 

K*(U10)-(K*°tt-)K+ 0.81 0.34 0.05 0.25 0.13 0.04 0.03 0.47 0.24 0.01 0.05 1.06 

K* a ~K*v S wave 0.83 0.22 0.20 0.12 0.27 0.06 0.06 0.64 0.01 0.01 0.05 1.13 

(t>p° S wave 1.34 0.35 0.41 0.29 0.05 0.10 0.10 0.22 0.23 0.02 0.23 1.53 

(j) P ° D wave 0.43 0.28 0.28 0.13 0.04 0.03 0.09 0.66 0.10 0.02 0.09 0.90 

(j){ir+TT-} s 0.43 0.57 0.19 0.16 0.09 0.10 0.16 0.19 0.00 0.05 0.19 0.83 

{K-n + } p {K+ir-} s 1.00 0.23 0.17 0.22 0.20 0.06 0.03 0.22 0.00 0.04 0.22 1.43 



region of the K s . The difference in parameter values between the nominal fit and that with 
the K S K + K~~ background fraction doubled is taken as the systematic uncertainty from this 
source. The largest background-related uncertainties are due to the statistical precision on 
the signal fraction. 

The uncertainty related to the modeling of experimental effects has three separate compo- 
nents: the acceptance, the resolution, and the mistag rate. The acceptance is incorporated 
in the fit using simulated data, as described in Sec. IIIIBl for CLEO III and CLEO-c data. 
To evaluate a systematic uncertainty related to the acceptance an alternative technique is 
used that is based on the product of individual particle efficiencies as a function of momen- 
tum in the laboratory frame. This is an almost identical procedure to the weighting used for 
the CLEO II.V normalization with simulated events. The only difference is an additional 
factor for CLEO-c data which is a function of tt + h~ invariant mass to account for the Kg 
veto. Such an approach is known to be simplistic compared to that using the fully simulated 
events given the integration over other variables on which the acceptance depends, such as 
polar angle, and the fact that it ignores correlations among the momenta. Therefore, it is 
considered a conservative approach to evaluating the systematic uncertainty related to the 
acceptance. The full difference in fit results for the two different techniques for incorporating 
the acceptance is taken as the systematic uncertainty. For the CLEO II.V data a uniform 
acceptance is assumed as an alternative model because the product of efficiencies is the 
technique used in the nominal fit. Most uncertainties due to the acceptance are around half 
a statistical standard deviation. 

The effects of resolution are ignored in the fit. The ensemble tests of simulated data 
reported in Sec. IIVCI are used to estimate the effect of resolution. The systematic variance 
due to the resolution is taken as the difference in the pull means squared for the fits performed 
with the generated and reconstructed four-momenta. The effect of resolution is found to be 
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very small or negligible on the fit parameters and fractions. 

The mistag rate for the different flavor-tag samples is varied within the uncertainties 
reported in Sec. [Til The resulting shift in the parameters is taken as the systematic uncer- 
tainty. 

The final source of systematic uncertainty considered is related to the overall performance 
of the fitter as demonstrated in Sec. IIV CI Small biases in the pull mean for some parameters 
are observed when fitting with the generated four-momenta. The largest bias is 40% of a 
statistical uncertainty with a significance of 3.5a from zero. The pull mean is assigned as 
a systematic uncertainty related to this fitter bias for each parameter. The source of such 
a bias could be the finite Monte Carlo statistics used to compute the normalization or the 
overall numerical precision of the fit. 

The systematic uncertainty is dominated by the mass and width assumed for the reso- 
nances; this leads to the measurements of the amplitudes and phases being systematically 
dominated. For the fit fractions the systematic and statistical uncertainties are approxi- 
mately of equal magnitude. 



E. CP violation search 

The baseline model is fit to the data allowing different values of Re(aj) and Im(aj) for 
D° and D° amplitudes. The fit is to the combined D° and D° flavor-tagged data to account 
correctly for the mistag rate. A comparison of the values of |a*| and fa for D° and D° 
amplitudes is given in Table IXII The systematic uncertainties are considered to be fully 
correlated between the D° and D° samples. The fit fractions are compared in Table Kill 
The fit fractions are used to determine the direct-CP asymmetry Acp, for each amplitude, 
which is defined as 

A = F d° ~ F W 
F D o + Fjjo ' 

where F^a and F-jyj are the fit fractions for D° and D° decays, respectively. There is no 
evidence for a significant CP asymmetry between any of the amplitudes. The sensitivity to 
Aqp varies between 5% to 30% among the amplitudes; therefore, the level of precision is 
not at the sub-percent level at which evidence of CP violation has been found in two-body 
D° decay @- 
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TABLE XI: |aj| and fa for the D° and -D amplitudes. The difference between the D° and D° 
parameter is also given in units of standard deviations (a). 



Parameter 


D° decays 


D° decays 


Difference (<r) 


#1(1270)- (K*°it-)K + \ ai \ 


0.35 ± 0.08 


0.39 ± 0.09 


0.3 


#i(1270)-(K* vr-)# + fa 


1.52 ±0.33 


0.98 ± 0.24 


1.3 


Ki(1270) + (p K + )K" \oi\ 


5.58 ±0.98 


5.96 ± 0.84 


0.3 


K x (1270)+ (p°K+)K- fa 


0.86 ±0.17 


0.71 ±0.15 


0.7 


K 1 (1270)-( P °K-)K + \oi\ 


7.03 ± 1.03 


6.39 ± 0.82 


0.5 


K 1 (127Q)-(p°K-)K+ fa 


0.41 ±0.19 


0.30 ±0.18 


0.4 


K*(UW) + (K*°tt + )K~ \oi\ 


5.39 ±0.91 


6.51 ±0.86 


0.9 


K*(U10) + (K* ir + )K~ 4>i 


0.99 ±0.18 


0.68 ±0.13 


1.4 


K*(U10)-(K*°ir-)K + \oi\ 


6.69 ±0.82 


6.75 ± 0.92 


0.0 


K*(UW)-(K*°ir-)K + fa 


1.38 ±0.20 


0.95 ±0.15 


1.7 


K*°K*° S wave \(n\ 


0.36 ± 0.05 


0.33 ± 0.04 


0.5 


K*°K*° S wave ^ 


0.47 ±0.17 


0.24 ±0.14 


1.1 


(j)p° S wave Oj 


1.02 ±0.16 


1.05 ±0.08 


0.2 


(j)p° S wave <j>i 


2.08 ±0.19 


1.77 ±0.07 


1.5 


4>p° D wave a« 


1.14 ±0.30 


1.69 ±0.23 


1.5 


D wave fa 


1.18 ±0.36 


1.48 ±0.14 


0.8 


0{vr+7r-} s 


5.76 ± 0.65 


6.22 ±0.65 


0.5 




1.75 ±0.19 


1.98 ±0.09 


1.1 


{K-TT+} p {K+TT~} s \0i\ 


84.7 ± 9.2 


81.7 ± 8.3 


0.2 


{K-TT+} p {K+7T-} s fa 


0.16 ±0.16 


0.13 ±0.10 


0.2 



TABLE XII: D° and D° fit fractions. The value of A C p is also given. 







Fit fraction (%) 


A C p 






D° Decays 


D° Decays 


(%) 


J FsTi(1270) + (#* ^ 


-)K- 


7.4 ± 1.1 


7.5 ±1.1 


-0.7 ±10.4 


#i(1270)-(#*°7T" 


-)K + 


0.9 ±0.4 


1.1 ±0.5 


-10.0 ±31.5 


K 1 (1270) + (p°K+)K- 


4.3 ± 1.1 


4.9 ±1.1 


-6.5 ± 16.9 


K 1 ll270)-(p°K )K + 


6.3 ± 1.1 


5.2 ± 1.0 


9.6 ± 12.9 


K*(U10) + (K*°tt- 


t)K- 


3.2 ±0.9 


4.8 ±1.0 


-20.0 ± 16.8 


#*(1410)-(#*°tt- 


~)K+ 


4.6 ±0.9 


4.7 ±0.9 


-1.1 ± 13.7 


K*°K*° S wave 




6.9 ± 1.2 


5.7 ± 1.2 


9.5 ± 13.5 


4>p S wave 




37.9 ±2.9 


40.0 ±2.9 


-2.7 ± 5.3 


4>p° D wave 




2.2 ±0.8 


4.8 ±1.2 


-37.1 ± 19.0 


(t>U + TT-} s 




9.0 ± 1.4 


10.7 ± 1.5 


-8.6 ±10.4 


{K-7T + } p {K+TT- 


} S 


11.3 ± 1.7 


10.7 ± 1.6 


2.7 ± 10.6 
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V. 7 SENSITIVITY STUDIES 



The first estimates of the sensitivity of B ± — > D (K + K^n + n^)K ± decays to the CP- 
violating parameter 7 were very promising [jjj]. These studies found that an uncertainty 
of approximately 10° is expected for a data set corresponding to an integrated luminosity 
of 2 fb _1 collected by LHCb. However, limited conclusions could be drawn as the model 
does not distinguish between D° and D° decays [H, @|. Therefore, these 7 sensitivity studies 
are repeated for the amplitude model presented in this paper, which is determined from 
flavor-tagged D° decays. 

An amplitude fit to simulated B — > D°(K + K~tt + tt~)K ± data is used to determine 
7. The fit is identical to that described in Ref. [3J and similar to the amplitude-model- 
dependent analyses of B + — > D(Kg7r + 7c~)K + j23[, which yield the most precise measure- 



ments of 7 to date |27l . |28| . The amplitude fit determines 7 from the distributions of 
D — > K + K~7t + tt~ events over four-body phase space, which are different for D mesons 
arising from B + or B~ decay. The distribution depends on the ratio between the CKM- 
suppressed and color-suppressed B~ — > D°K + amplitude and the Cabibbo-favored and 
color-allowed B — > D°K~ amplitude and is parametrized as rse^ 5s_7 \ where r# is the 
magnitude of the amplitude ratio and 5 b is the CP-invariant strong-phase difference between 
the amplitudes. The values of rs and 5b are also determined by the amplitude fit. 

An ensemble study of 200 simulated data sets containing 2000 B ± — > D(K + K~'k + 'k~)K ± 
events each, split evenly between the B meson charges, is used to estimate the sensitivity 
to 7. The number of events in each data set corresponds approximately to that expected in 



a few years running of LHCb 29]. The values of 7, r^, and 5b are assumed to be 70°, 0.1 



and 130°, respectively, which are approximately the world average values [30j. These data 
samples are then fit assuming the amplitude model used in the generation, with 7, t\b, and 
5b as free parameters. The average uncertainty on 7 from the ensemble of experiments is 
(11.3 ± 0.3)° and the pull distribution formed from the fitted and generated values follows 
a normal distribution, indicating the results are unbiased. The pull distributions for r# 
and 5b are also normal. Ensemble studies for the alternative models given in Appendix O 
are also performed. Most alternative models yield a sensitivity to 7 similar to that of the 
baseline model. However, model 6 leads to an average uncertainty approximately 50% worse 
than the other models. Similar model-dependent variations in uncertainty are reported in 
Ref. Q. However, overall these results confirm that there is significant sensitivity to 7 in 
B ± b(K + K-TX + -n-)K ± decays. 



VI. CONCLUSIONS 

The first amplitude model for D° — > K + K~7t + tt~ decay derived from flavor-tagged data 
has been presented. The data used are from e + e~ collisions at center-of-mass energies close 
to cc threshold and in the region of the T resonances. CP-tagged quantum-correlated data 
recorded at the -^(3770) resonance are also used in the fit. The model indicates that the 
quasi-two-body decay D° — > <pp° is dominant, with significant contributions from the fol- 
lowing intermediate states: D° ^(1270)^^, D° K^UIO^K^, and D° <j)n + n-. 
In addition, there is a significant D° — > K*°K*° contribution indicating that final-state 
interactions play a significant role in the decay. There is also a non-resonant contribution 
of around 10%; the best fit to data is achieved when there is relative angular momentum 
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TABLE XIII: Spin Factors (S) for various amplitudes. In the decay chains, S, P, V and A stand 
for scalar, pseudoscalar, vector, and axial vector, respectively. Letters in square brackets refer to 
whether the decay products are in a relative S, P, or D wave state of relative orbital angular 
momentum. Spin factor number [3] with D[D] — y Vi V® actually corresponds to a superposition of 



D and S wave, following the choice of basis in Ref. 31]. If no angular momentum is specified, 
the lowest angular momentum state compatible with angular momentum conservation and, where 
appropriate, parity conservation, is used. 

Number Decay chain S 

1 D APi, A[5] => VP 2 ,V => P3P4 pfpjA^Z^Vf 

2 d[ ( s]^v 1 v 2 ,v 1 ^p 1 p 2 ,v 2 ^p 3 p 4 z 1 (y 1 ) 1/ z 1 (y 2 y 

3 D[P]^ViV 2 ,Vi^PiP 2 ,V 2 ^P 3 P 4 taMSPD&lv^k 

4 B[D] -> V 1 Y 2 ,Y 1 -> P X P 2 , V 2 -> P3P4 Z^Vija pfc < 

5 D^VS,V^PiP 2 ,S^P 3 P 4 Ps z i( v )n 

6 D^V 1 P 1 ,Vi^V 1 P 2 ,V 2 ^P 3 P4 ^M&PvAyp^k 



among the particles in this contribution. The accuracy of the model parameters is limited 
by the uncertainties on the i^i(1270) _ resonance parameters. 

The amplitude model presented has been used to search for CP violation in the decay 
by determining the fit fractions separately for D° or D° decays. The fit fractions are found 
to agree within the uncertainties, indicating no CP violation in the decay at the level of 
a few percent. The amplitude model has also been used in a sensitivity study of B ± — y 
D°(K + K~7t + 7t~)K ± decays to the CP-violating parameter 7. The study indicates that 7 
can be determined with a precision of (11.3 ± 0.3)°, assuming the baseline model, using this 
decay at LHCb. A similar precision would be expected at future flavor facilities. 



Acknowledgments 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Natural Sciences and Engineering Research 
Council of Canada, and the U.K. Science and Technology Facilities Council. 



Appendix A: Spin factor definitions 

The spin factors used are those calculated in [31]. These are expressed in terms of the 
four-momenta of the particles involved. For the decay R — y A, B the following notation is 
used: 

• Pr,Pa, and vb are the four-momenta of the resonance R, its decay products A and, 
B, respectively. This are related such that 

Pr = Pa+Pb- (Al) 

• qn represents the difference of the four-momenta between the decay products of R: 

Qr = Pa - Pb ■ (A2) 
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Here the ordering of the particles is important; the momentum of the second daughter 
particle listed in the decay chain is always subtracted from the momentum of the first. 



For the purpose of representing the spin factors in a concise way, the following functions 
are defined: 



P(p, m y« = g ^-UL and 

mr (A3) 



Z l (q,p,m 2 ) a = qtl P(p,^ a 



m 



For operators related to the resonance R the expressions are simplified by using the defini- 
tions 



P(RY a = P(p R , m 2 R Y a and 
Z x {R) a ee Zi(q R ,p R ,m 2 R 



(A4) 



For non-resonant contributions m R is replaced with p 2 , as in P{p,p 2 ) tia for example. The 
spin factors used are listed in Table IXIII1 

It is clear that the exact matching of the particles P±, P2, P3, and P\ in the spin factor 
definition to the final state particles in the decay is important, as many spin factors change 
sign under swapping a pair of particles due to terms such as q R = p\ —p 2 - For the amplitudes 
used in this paper, the particle ordering is given in Table IXIVI 



Appendix B: Amplitudes tested 

Below is a list of all the different amplitudes that are tested when determining the best 
model. For final states that are flavor specific the charge-conjugate amplitude is not listed 
but is one of the amplitudes that is tested. 

• Cascade amplitudes containing a higher K* resonance 

- fs: 1 (1270) + (ir* 7r+)fs:-, K x (1 270) +(K*( 1430) it + )K~, ATi(1270) + (p° 'K+)K~ , and 
Ki(1270) + (uK + )K- 

- ir 1 (1400) + (K*°7r + )ir- 

- K 2 {IA^) + {K*°tx + )K- and K 2 (U30) + (p°K+)K- 

- K*(l680) + (K* ir + )K~ and K*(lQ80) + (p°K+)K- 

• Quasi-two-body amplitudes 

- K*°W° S, P, and D wave 

- 0p° S, P, and D wave 

- (f)0J S wave 

- 0/2(127O) P and D wave 

• Single resonance amplitudes 

- p{K + K~} s ; p{K+K-} P S, P, and D wave; and p{K + R-} D P and D wave 
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TABLE XIV: Spin factors used for different decay chains, including the particle numbering scheme. 
The 2 nd column refers to the spin factors as numbered in Tab IXIlH and the particles Pi, P2, P3, 
and P4 refer to the numbers as defined in Tab IXIIII 



Decay chain 


Spin factor 


number Pi 


P2 


P3 


P4 


ifi(1270) + (if * u 7r + )if 


rn 

1 


if- 


_|_ 


if + 




ifi(1270)~ (if*°7r~)if+ 


■ = ■ 

1 


if+ 


7T~ 


if" 


7T+ 


i^i(1270) + (pK + )K~ 


1 


if~ 


K + 


_(_ 


7T~ 


K 1 (1270) (pK )if + 


1 


K ^ 


T ' 

K 




4. 


if*(1410) + (if* u 7r+)if 





if- 


1 

7T+ 


K + 


7T~ 


if* (1410) ~ (K*°ir~)K + 


6 


K + 


7T~ 


K~ 


4. 
7T + 


if * u if *° 5 wave 


2 


K + 


7T~ 


K~ 


4- 

7T + 


if *°if *° P wave 


3 


if+ 


7T~ 


K~ 


7T+ 


if* {if-yr+jp P wave 


i 


if+ 


7T~ 


K~ 


7T+ 


if*°{K + 7r-} p wave 


3 


if+ 


7T~ 


K~ 


7T+ 


if *°if *° D wave 




if+ 


7T~ 


if- 


7T+ 


{A'-7r+} p {K+7r-}p D wave 


i 


K+ 


7T _ 


if- 


7T+ 


4>p° S wave 


2 


if+ 


if" 


7T+ 


7T~ 


0{7r + 7r _ }p S wave 


= 

2 


if+ 


if~ 


1 

7T+ 


7T~ 


p {if+if~}p 5 wave 


2 


K+ 


if- 


7T+ 


7T~ 


(f){ir + ir~} p P wave 


3 


if+ 


if" 


7T+ 


7T~ 


0/?° D wave 


i 


if+ 


if- 


7T+ 


7T~ 


</>{7r + 7r~}p P wave 


a 


if+ 


K~ 


7T+ 


7T~ 


{ir + ir-} p {K + K-} p D wave 


si 


if+ 


K~ 


7T+ 


7T~ 


(/>{7r+,vr-} 5 


5 


if+ 


if- 


7T+ 


7T~ 


{ET-7r + }p{K+^-} s 


5 


K~ 


vr+ 


if+ 


7T _ 


{K+7r-} p {K-7r+} s 


El 


if+ 


7T- 


if- 


7T+ 



- if*°{if-7r + } s ; if*°{if-7r+}p S, P, and P> wave; and K *°{K-n + } D P and D 
wave 

- 0{7T + 7r-}5; (/){7r + 7r-}p S, P, and .D wave; and 0{7r + 7r-}p, P and P> wave 

- io(980)°{7r+7r-} 5 and / (980)°{if+if-} s 

- / 2 (1270)°{if+if-} s 

- u{K+K-} s 

• Non-resonant amplitudes 

- {if+if-} s {7r+7r-} s ; {K + K-} s {7i+7i-}p; {if +if-} P {7r+7r-} 5 ; 
{if + if-}p{7r+7r-}p S, P, and D wave; {if + if-} 5 {7r+7r-} D ; 
{if + if-}p{7r+7r-} 5 ; {if+if-}p{7r+7r-} D P and P> wave; {if + if-}p{7r+7r-}p 
P and D wave 

- {if + 7r-} 5 {if-7r+}p; {A^+7r-}p{if-7r+} 5 ; {if+7r-}p{if-7r+}p S, P, and D 
wave; {K + 7i-} s {K-tt + } d ; {K + 7t-} d {K-tt + } s ; {if+vr-}p{if-7r + } D P and D 
wave; {if + 7r-}p{if-7r + }p P and D wave 
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TABLE XV: Fit fractions (in %) and \ 2 l u f° r alternative models in a combined fit to all flavor and 
CP-tagged data sets. The uncertainties are statistical. Where it is necessary to specify the angular 
momentum state of pairs of non-resonant particles, this information is given in the subscript. The 
baseline model adopted in this paper is number 1. 



Model 


1 


2 


3 


4 


5 


6 


7 


8 


9 



K 1 (1270)+(K* u tt+)K- 7.3±0.8 7.0±0.8 7.0±0.8 7.0±0.8 7.0±0.8 8.7±0.9 6.7±0.8 8.4 ±1.0 8.7 ± 0.9 

#1(1270)- (K*°tt-)K+ 0.9±0.3 0.7±0.3 0.8±0.3 0.8±0.3 0.9±0.3 1.7±0.5 0.8±0.3 1.7 ±0.5 0.8 ± 0.4 

Ki(1270) + (p°K+)K- 4.7±0.7 4.7±0.7 4.8±0.7 4.6±0.7 6.0±0.8 5.1±0.7 6.8±1.1 6.0 ± 0.9 4.5 ± 0.8 

Ki(1270)- (p°K-)K+ 6.0±0.8 5.8±0.7 5.9±0.7 6.0±0.8 3.7±0.8 6.3±0.8 3.5±0.8 7.9 ±1.3 5.9 ± 0.9 

K* (1410)+ (K*°tt+)K- 4.2±0.7 4.2±0.7 3.4±0.7 4.2±0.7 4.2±0.7 4.1±0.7 3.2 ± 0.7 3.2 ±0.7 

K* (1410) " (K* %- ) K+ 4.7±0.6 4.5±0.6 2.9±0.6 4.5±0.6 4.7±0.6 4.5±0.6 5.1 ±0.8 4.7 ± 0.7 

K*° {K-tt+} p P wave 3.9±0.7 — 

K*° {K+ir~ } p P wave 3.0±0.5 

K*°~K*° S wave 6.1±0.8 6.1±0.8 6.1±0.8 6.1±0.8 6.3±0.8 4.2±0.6 6.3±0.8 3.6 ± 0.8 4.4 ± 0.8 
K*°~K*0 p wave 

^{7r+7r-} s 10.3±1.0 9.0±0.9 9.1±0.9 9.6±1.0 10.1±0.9 8.9±0.9 7.3 ± 0.9 9.1 ±0.09 

4> {w+w- } p S wave 6.3±1.1 6.1±1.1 11.4±2.5 5.9±1.1 

4>{tt+tt-} p P wave — 4.7±0.6 — 



6.1±0.8 


6.1±0.8 


6.1±0.8 


6.1±0.8 






0.9±0.3 




10.3±1.0 


9.0±0.9 


9.1±0.9 


9.6±1.0 




6.3±1.1 


6.1±1.1 


11.4±2.5 








2.0±1.1 


38.3±2.2 


21.0±2.0 


21.2±1.9 


15.9±2.6 


3.4±0.7 


3.7±0.7 


3.7±0.7 


1.1±0.6 


10.9±1.2 


10.9±1.1 


11.0±1.1 


10.8±1.1 



4> {77+77 }p D wave 

</>p° S wave 38.3±2.2 21.0±2.0 21.2±1.9 15.9±2.6 38.0±2.0 26.8±1.3 20.9±1.9 34.6 ±2.1 36.7 ±0.2 

4>p° D wave 3.4±0.7 3.7±0.7 3.7±0.7 1.1±0.6 3.5±0.7 3.8±0.7 2.9 ±0.7 2.8 ± 0.7 

p° {K+K~} p S wave 1.2±0.6 1.2±0.6 

{K-tt+} p {K+tt-} s 10.9±1.2 10.9±1.1 11.0±1.1 10.8±1.1 9.3±1.1 9.1±1.1 9.0±1.2 

{K+tt-} p {K~n+} s 11.6 ±1.5 5.0 ±1.0 

{K-tt+} p {K+tt~} p D wave 13.8±1.2 

{7r+7T-} p {K+K~} p D wave — — — — — 2.6±0.7 — — — 

Sum 96.7±2.6 84.0±2.0 82.9±2.0 84.0±2.0 94.9±2.3 80.9±1.4 82.5±2.0 92.4 ± 2.4 94.9 ± 2.4 

X z /v 1.63 1.64 1.65 1.65 1.66 1.66 1.66 1.77 1.78 



Appendix C: Alternative models 

The fit fractions and x 2 l v f° r tb_ e best seven fits to the D° — > K + K~tt + 7t~ (models 1 
to 7) data are shown in Table IXVI Model 1 is chosen as baseline. Models 8 and 9 contain 
different variations of the non- resonant amplitude component. 



Appendix D: Accounting for quantum correlations 

In the analysis the CLEO-c flavor-tagged data are assumed, apart from the proportion 
of events that were mistagged, to form a pure flavor sample. However, since the D°D° pair 
results from the decay of the J PC = 1 ^(3770) particle, the D mesons are produced in 
a correlated state. For example, when both the D° and D° decay to CP-eigenstates, these 
states will have opposite CP. The antisymmetric wave function which describes the decay 
is 

[AooisjA, - A 3 A w (s)] . (Dl) 

Here, Aj (Aj) is the amplitude for the decay of the other D° (D°) in the decay. In this 
paper the CLEO-c flavor tagging is provided by tagging the non-signal D in its inclusive 
decay to kaons under the assumption of a Cabibbo-favored decay. For example D° — > K~X 



-ttO 



and D — > K + X. Therefore it is useful to consider the case where state j of Eq. fIDll) is 
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denned to be a specific tag, labeled by the subscript t, so that 

Aj = A(D° -> K~X t ) = K t and (D2) 

Aj = A(D° K~ = K t ^e i5< , (D3) 

where K t , k t , and 5t are real numbers. Therefore, from Eq. ( 1D1[) the full decay rate for the 
decay of interest against the specific tag defined by Eq. (1D3I) . will be given by 

r t oc K*(\A D o\ 2 + k 2 t \Aw\ 2 

- 2k t [Re[A D oA- m ] cos S t + lm[A D oA w ] sin 5 t }) . (D4) 

summing over all possible tags t, which will all in general have different values for k t , K t , 
and S t . Hence, the full decay rate will be given by 

r oc (^ 2 )[l^| 2 +(^gi) \a w \ 2 

r , 11 f^tK?k t cos S t \ r i in f^tKfkt sin &\ , . _ , 

- 2(Re[^o^] ^ ^ « j + Im[^o^] ^ M )] . (D5) 

Defining 

W - (^f ) , 

// jt\ - f^tKfk t cos 5 t \ 

(k cos 5) = I ^-^2 1 , and (D7) 

(A; sin 5) = I 1 , (D8) 

Eq. ( 1D5I) can be rewritten as 

T oc (Et^O^ol' + I^I^A; 2 ) (D9) 

- 2 (Re[^ D o^^o](A;cos5) + Im[i D o^] (k sin 5))]. 

Since the tags in question are dominated by Cabibbo-favored and doubly-Cabibbo- 
suppressed decays it is noted that k t « 0.05. Therefore, the term \A-po\ 2 (k 2 ) can be ne- 
glected. It is also assumed that the parameters relating to the tags - (J2 t K 2 ) (kcos6), and 
(A; sin S) - will not vary over the K + K^-k + -k^ Dalitz space; this assumption allows the term 
(E t K 2 ) to be absorbed into the normalization, leaving only the parameters (kcos5) and 
(A; sin 5) to be determined from the data. 

A fit to CLEO-c 3770 flavor-tagged data is performed with the distribution described by 
Eq. (D9), with (kcos5) and (k sin 5) as additional free parameters. The results of this fit 
are (A; cos 5) = 0.061 ± 0.042 and (k sin 5) = 0.029 ± 0.007. The difference in the values of 
Re(aj) and Im(aj) from the quantum-correlated fit compared to the nominal fit are used to 
estimate the systematic uncertainty. 

Quantum correlations are also present for data produced with a center-of-mass energy 
of 4170 MeV. However the D° and D particles are generally not produced directly at this 
energy but via the decay of higher mass resonances; therefore, the interference effects are 
expected to not be the same as for the CLEO-c 3770 data. The measured values of (k cos 5) 
and (k sin 5) are 0.051 ± 0.032 and 0.037 ± 0.011, respectively. 

The systematic shifts due to quantum correlations at the ^(3770) and -0(4170) center- 
of-mass energies are added in quadrature in order to give the total systematic shift due to 
these effects. 
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Appendix E: EPAPS information for the paper "Amplitude Analysis of D —> 

K+K--K+1T-" 

1. Additional information about the baseline fit 

Table IXVII gives a key to the amplitude and parameter indexing used in this addendum. 
Table IXVIII gives the statistical correlation matrix among the fit parameters for the best-fit 
model (model 1) described in the main body of the text. Table IXVlHI gives the interference 
fractions for best-fit model. The interference fraction between the m th and n th amplitudes, 
Z m>n , is defined as: 

_ j 2Re(A m A* n )d Pj 

J |*4 D o| 2 <% 

The interference fractions are listed in order of their absolute magnitude and only terms 
with a magnitude greater than 0.5% are listed. 



TABLE XVI: Amplitude and parameter indexing used in this addendum. 



Amplitude 




Index 


J PsTi(1270) + (K* 7r- 




1 


i^i(1270)-(K*°7r- 


-)K+ 


2 


K 1 {mQ) + {p Q K + 


)K- 


3 


Ki{l27Q)-{p°K- 


)K + 


4 


K*(U10) + (K*°tt- 


^)K- 


5 


K*(U10)-(K*°tt- 


-)K+ 


6 


K*°K*° S wave 




7 


4>p° S wave 




8 


4>p° D wave 




9 


(p{TT + TT-} s 




10 


{K-TT + } p {K + 7T~ 




11 



2. Fit parameters for alternative models 

Tables IXIXI and IXXI contain the real and imaginary parts of fitted for alternative 
models 2-5 and 6-9, respectively. 
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Im(<22) Re(ci3) Im(a3) Re(a4) Im(a4) Re(as) 1111(05) Rc(ag) Im(ag) Rc(a7) 1111(07) Rc(as) Im(ag) Rc(ag) Im(ag) Rc(aio) Im(aio) Rc(an) Im(an) 



Rc(a 2 ) 


0.164 


-0.225 


0.361 


0.022 


0.089 


0.155 


0.214 


0.455 


0.048 


-0.212 


0.408 


-0.347 


-0.120 


0.236 


0.042 


0.297 


-0.246 


0.002 


0.409 


Im(ct2) 


1.000 


-0.223 


0.070 


-0.218 


0.092 


0.019 


0.280 


0.196 


0.176 


-0.334 


0.080 


-0.036 


-0.177 


0.038 


0.175 


0.043 


0.188 


-0.175 


0.188 


Re (013) 




1.000 


-0.409 


-0.057 


-0.139 


-0.016 


-0.475 


-0.339 


-0.367 


0.380 


-0.362 


0.359 


0.458 


-0.277 


-0.272 


-0.441 


-0.080 


0.423 


-0.460 


Im(a3) 






1.000 


0.269 


0.076 


0.318 


0.167 


0.542 


0.008 


0.139 


0.520 


-0.335 


0.078 


0.252 


-0.164 


0.308 


-0.420 


0.044 


0.534 


Re (04) 








1.000 


-0.261 


0.204 


-0.260 


0.113 


-0.336 


0.444 


0.079 


-0.178 


0.203 


0.094 


-0.252 


0.085 


-0.474 


0.272 


0.043 


Im(a4) 










1.000 


0.220 


0.188 


0.290 


0.099 


-0.061 


0.398 


0.072 


0.180 


-0.120 


-0.015 


-0.064 


-0.011 


0.206 


0.284 


Re (05) 












1.000 


0.024 


0.215 


-0.039 


0.189 


0.304 


-0.107 


0.217 


0.025 


-0.191 


0.059 


-0.307 
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Im(as) 
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0.368 


0.353 


-0.472 


0.401 


-0.367 


-0.526 
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Re(a 6 ) 
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-0.194 
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0.084 
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Im(a 6 ) 
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0.216 


0.340 
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0.554 


-0.321 
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Im(a9) 






























1.000 


0.224 


0.302 


-0.381 


0.144 



TABLE XVIII: Interference fractions for the best-fit model in %. Only interference fractions greater 
than 0.5% are given. 





Value (%) 
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2^4,8 


5.7 


Zij 


3.3 


2^3,8 


3.0 


2^3,4 


2.8 


2^10,11 


-2.6 


2^1,3 


1.5 


2^3,7 


1.3 


2^7,8 


1.1 


2^5,6 


1.0 


2^1,8 


1.0 


2^4,9 


-0.9 


2^4,7 


0.8 


2^3,9 


-0.6 
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